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ABSTRACT 

We present a model for the distribution of void sizes and its evolution in the context of 
hierarchical scenarios of gravitational structure formation. We find that at any cosmic 
epoch the voids have a size distribution which is well-peaked about a characteristic void 
size which evolves self-similarly in time. This is in distinct contrast to the distribution 
of virialized halo masses which does not have a small-scale cut-off. 

In our model, the fate of voids is ruled by two processes. The first process affects 
those voids which are embedded in larger underdense regions: the evolution is effec- 
tively one in which a larger void is made up by the mergers of smaller voids, and is 
analogous to how massive clusters form from the mergers of less massive progenitors. 
The second process is unique to voids, and occurs to voids which happen to be em- 
bedded within a larger scale overdensity: these voids get squeezed out of existence as 
the overdensity collapses around them. It is this second process which produces the 
cut-off at small scales. 

In the excursion set formulation of cluster abundance and evolution, solution of 
the cloud- in- cloud problem, i.e., counting as clusters only those objects which are 
not embedded in larger clusters, requires study of random walks crossing one-barrier. 
We show that a similar formulation of void evolution requires study of a two-barrier 
problem: one barrier is required to account for voids-in-voids, and the other for voids- 
in-clouds. Thus, in our model, the void size distribution is a function of two parameters, 
one of which reflects the dynamics of void formation, and the other the formation of 
collapsed objects. 

Key words: galaxies: clustering - cosmology: theory - dark matter. 



1 INTRODUCTION 

An overwhelming body of observational and theoretical ev- 
idence favours the view that structure in the Universe has 
risen out of a nearly homogeneous and featureless primor- 
dial cosmos through the process of gravitational instability. 
Almost all viable existing theories for structure formation 
within the context of this framework are hierarchical: the 
matter distribution evolves through a sequence of ever larger 
structures. 

Hierarchical scenarios of structure formation have been 
successful in explaining the formation histories of gravita- 
tionally bound virialized haloes. They provide a basic frame- 
work within which more intricate aspects of the formation 
of a wide range of cosmic objects, ranging from galaxies to 
rich clusters, may be investigated. In particular, a fully an- 
alytical description of the collapse and virialization of over- 
dense dark matter halos has been developed. The approach, 
originally proposed by Press & Schechter (1974), and later 
modified by Epstein (1983) and Bond et al. (1991), has led 



to simple and accurate models for the abundance of massive 
haloes which results from hierarchical gravitational cluster- 
ing. This framework has come to be called the excursion set 
approach. 

The excursion set approach provides a useful framework 
for thinking about the formation histories of gravitationally 
bound virialized haloes in scenarios of hierarchical structure 
formation. It provides analytic approximations for the dis- 
tribution of halo masses, merger rates, and formation times 
which are quite accurate (Lacey & Cole 1993), and can be 
extended to provide estimates of the distribution of the mass 
in randomly placed cells (Sheth 1998). A key ingredient in 
the original approach, inherited from the pioneering work 
of Press & Schechter (1974), is the assumption that virial- 
ized objects form from a smooth spherical collapse. In real- 
ity the collapse can be quite different from spherical; recent 
work has shown that ellipsoidal collapse can be incorporated 
into the approach, with reasonable improvements in accu- 
racy (e.g. Sheth, Mo & Tormen 2001). 

Models based on spherical evolution are difficult to rec- 
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oncile with the spatial patterns which characterize the cos- 
mic matter distribution. The observed world of galaxy red- 
shift surveys, and the artificial world of numerical simula- 
tions of cosmic structure formation, are both characterized 
by filamentary and sheetlike structures. Such weblike pat- 
terns represent distinctly non-virialized structures for which 
gravitational contraction of initially aspherical density peaks 
has only been accomplished along one or two dimensions. At 
first sight, such weblike configurations would seem to be be- 
yond the realm of the idealized excursion-set description. 

Nevertheless, in this study we show that the formation 
and evolution of foamlike patterns can indeed be described 
by the excursion set analysis. This is accomplished by focus- 
ing on the evolution of underdense regions, the voids, rather 
than overdensities in the matter distribution. Whereas much 
of the mass in the universe is bound-up in virialized struc- 
tures, most of the volume is occupied by large underdense 
voids: voids are the dominant component of the Megaparsec- 
scale galaxy and matter distributions. In a void-based de- 
scription of structure formation, matter is squeezed in be- 
tween expanding voids, and sheets and filaments form at 
the intersections of the void walls (Icke 1984; Van de Wey- 
gaert 1991, 2002). Such a view is supported by Regos & 
Geller (1991), Dubinski et al. (1993) and Van de Weygaert 
& Van Kampen (1993), who give clear and lucid descrip- 
tions of how voids evolve in numerical simulations of gravi- 
tational clustering. We will stick to this basic framework in 
the present study. 

We will argue that low density regions are the objects- 
of-choicc for working out a succesful analytical description 
of cosmic spatial structure, if it is to be based upon the ide- 
alization of spherical symmetry. This is because, in many 
respects, voids are ideally suited for an excursion set analy- 
sis based on a spherical evolution model. This is despite the 
fact that voids form from negative density perturbations in 
the initial fluctuation field, and neither maxima nor minima 
in the primordial Gaussian field, are spherical (see Bardeen 
ct al. 1986). However, in marked contrast to the evolution 
of density peaks, primordial asphericity of negative density 
perturbations is quickly lost as they expand: the generic evo- 
lution is towards an approximately spherical tophat geome- 
try (Icke 1984). Moreover, the velocity structure of uniform 
density voids is simple to understand; an observer in the in- 
terior will observe a Hubble-type velocity field. All of this 
is discussed in some detail in Appendix A, which describes 
the evolution of a single isolated void. 

Although the image of a large scale matter distribu- 
tion organized by expanding voids is appealing, in its basic 
form, the description essentially involves an extrapolation 
of single void characteristics to an entire random popula- 
tion of strictly distinct and non-interacting peers, each of 
them undisturbed smoothly expanding bubbles. This dis- 
cards one of the most crucial and characteristic aspects 
of cosmic structure formation — that there are no isolated 
voids, nor smoothly unstructured ones. Any complete anal- 
ysis will have to take into account the complications which 
arise from 

• the substructure present within the primordial volume 
occupied by the void, and 

• the inhomogeneous matter distribution in its vicinity. 

The existence of internal void structure is not unexpected. 



The void shown in Figure 1, selected from a large N-body 
simulation of cosmic structure formation, shows the exis- 
tence of structure on all scales. The figure shows three suc- 
cessive zoom-ins on the inner parts of the void; all exhibit 
some measure of internal structure, although substructure 
is less pronounced in the emptiest inner regions. 

As was mentioned above, all viable cosmological struc- 
ture formation scenarios imply a hierarchical mode of struc- 
tural growth. The formation of any object involves the fusion 
of all substructure present within its realm, including the 
small-scale objects which had condensed out at an earlier 
stage. Underdensities are organized similarly — in the evo- 
lution of a void we may identify two, intimately related, 
processes: 

• a bottom-up assembly, in which a void emerges as a ma- 
ture and well-defined entity through the fusion and gradual 
erasure of its internal substructure, and 

• the interaction of the void with its surroundings, mark- 
ing its participation in the continuing process of hierarchical 
structure formation. 

Considerable insight into the evolution of voids came from 
the rigorous and insightful study by Dubinski et al. (1993). 
Following an analytical study of (isolated) spherically sym- 
metric voids by Blumenthal et al. (1992), they used N-body 
simulations to study the evolution of the void hierarchy from 
a set of artificial and simplified initial conditions, consist- 
ing of various levels of hierarchically embedded spherical 
tophat voids. They showed that adjacent voids collide, pro- 
ducing thin walls and filaments as the matter between them 
is squeezed. Mainly confined to tangential motions, the pe- 
culiar velocities perpendicular to the void walls are mostly 
suppressed. The subsequent merging of voids is marked by 
the gradual fading of these structures while matter evac- 
uates along the walls and filaments towards the enclosing 
boundary of the "void merger" . The timescale on which the 
internal substructure of a void is erased is approximately 
the same as that when the void itself approached "nonlin- 
earity" (Appendix A gives a precise definition of what is 
meant by nonlinearity) . At nonlinearity, smaller-scale voids 
collide and merge with one another, effectively dissolving 
their separate entities into one larger encompassing void. 
Only a faint and gradually fading imprint of their original 
outline remains as a reminder of the initial internal substruc- 
ture. As this (re)arrangement of structure progresses to ever 
larger scales, the same basic processes repeat. 

N-body simulations of voids evolving in more generic 
cosmological circumstances by Van de Weygaert & Van 
Kampen (1993) (also see Van de Weygaert 1991) yielded 
similar results. This prompted them to suggest the existence 
of a natural void hierarchy, in which small-scale voids em- 
bedded within a pronounced large-scale void gradually fade 
away. An illustration of such a void hierarchy process, within 
the context of the CDM scenario, is shown in Figure 2. The 
major characteristics of an evolving void hierarchy, the grad- 
ual blending of small-scale voids and structures into a larger 
surrounding underdensity, is clearly visible in the sequence 
of six timesteps. 

However, the artificial arrangement of voids embedded 
within voids represents only one aspect of reality — it misses 
a crucial component of the development of a void hierarchy. 
An evolving void hierarchy not only involves the merging of 
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Figure 1. Spatial structure in a void-like region selected from an N-body simulation of structure formation in the SCDM scenario. 
Three consecutive zoom-ins centered on the core of the void are shown: a 45Mpc diameter particle sphere (top), 36Mpc (bottom left) 
and 30Mpc (bottom right). The existence of substructure within the void region is readily apparent, although it becomes more faint and 
tenuous towards the increasingly depleted interior of the void. 



small voids into larger voids, but also the disappearance of 
small voids as they become embedded in larger-scale overden- 
sities. Thus, in contrast to the process of dark halo forma- 
tion, the emerging void hierarchy is ruled by two processes 
instead of one. The main goal of this paper is to incorpo- 
rate both processes into a model of the void hierarchy. We 
do this by combining the spherical evolution model with the 
excursion set approach. When used to describe the evolu- 
tion of overdense clouds, the excursion approach requires 
consideration of a one-barrier problem, the single barrier 
representing what is required for collapse in the spherical 
evolution model. We show that the excursion set formulation 
of the void hierarchy requires consideration of two-barriers: 
one barrier is associated with the collapse of clouds, and the 
other with the formation of voids. The resulting framework 
is able to describe realistic settings of random density fields 
in which voids interact with their surroundings. 

This paper is organized as follows. Section 2 discusses 
important generic properties of isolated voids, which grow 
from depressions in the primordial density field, propelled 
by the perturbed gravitational field. The spherical model 



forms the core of further analytical considerations, and is 
discussed in some detail in Appendix A. Section 3 discusses 
the generic effects of larger scale stucture on the evolution of 
voids. Two crucial processes which shape the void hierarchy 
are described: the void-in-void mode and the void-in-cloud 
mode. How these processes can be incorporated into the ex- 
cursion set approach using two-barriers is the subject of Sec- 
tion 4. Section 5 describes the associated distribution of void 
sizes, which is predicted to have a universal form, and to be 
peaked around a characteristic value. One of the results of 
Section 5 is to show that peak-based models should be rea- 
sonably accurate for the largest voids, but, because they 
account neither for the void-in-void mode nor for the void- 
in-cloud mode, they predict many more small voids than 
does the excursion set approach. Appendix B discusses the 
"basic troughs model", which assumes there is a one-to-one 
identification between minima in the primordial Gaussian 
density field, with centres of voids in the evolved (and non- 
linear) matter distribution. This also serves to define nota- 
tion for the "adaptive troughs model" which is described in 
Section 5.1. 
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Figure 2. Void evolution. Six timesteps in the evolution of a void region in a 128 3 particle N-body simulation of structure formation in 
an SCDM model: top-left to bottom-right shows expansion factors a cxp =0.1, 0.2, 0.3, 0.35, 0.4 and 0.5 (the present time has a oxp = 1.0). 
Initial conditions were defined such that they would focus in on a 3<r(4ft- 1 Mpc) void, using a constrained random field code (Van de 
Weygaert & Bertschinger 1996). The sequence shows the gradual development of a large void of diameter xs 25h~ -'-Mpc as the complex 
pattern of smaller voids and structures which had emerged within it at an earlier time, merge with one another. This illustrates one 
aspect of the evolving void hierarchy: the void-in-void process. 



Section 6 presents various other aspects of the hierar- 
chically evolving void population. Global parameters, such 
as the fraction of mass in the cosmos contained within void 
regions, along with the fraction of space occupied by voids, 
are readily derived from the void size distribution. In addi- 
tion, the formalism is applied towards a reconstruction of 
the ancestral history of a given void, followed by an evalu- 
ation of the environmental influence on basic void proper- 
ties. We also put forward suggestions towards an analytical 
treatment of the influence of the void environment on the 
galaxies that may form within. Finally, we indicate how an 
assessment of the evolution of dark matter clustering may 
be predicated on our formalism. In Section 7, we provide 
an overview of our results and seek to embed these in the 
wider context of the study of hierarchical structure forma- 
tion. We also comment on how our results for the distribu- 
tion of voids in the dark matter distribution may be related 
to observations of voids in the galaxy distribution. Although 
our model provides a useful framework, developing a more 
detailed model is beyond the scope of this work. The results 
of numerical studies of void galaxies in semi-analytic galaxy 
formation models are described by Mathis & White (2002) 
and Benson et al. (2003). 



2 EVOLUTION OF ISOLATED VOIDS 

The basic features of voids can be understood in terms of 
the evolution of isolated density depressions. The net density 
deficit brings about a sign reversal of the effective gravita- 
tional force: a void form from a region which induces an 
effective repulsive peculiar gravity. 

In physical coordinates, overdense regions expand 
slightly less rapidly than the background, reach a maximum 
size, and then turn around and finally collapse to vanishingly 
small size (this is strictly true only in an Einstein de-Sitter 
or closed Universe). In contrast, underdense regions will not 
turn around: they undergo simple expansion until matter 
from their interior overtakes the initially outer shells. The 
generic characteristics of these evolutionary paths may be 
best appreciated in terms of the evolution of isolated spher- 
ically symmetric density perturbations, either overdense or 
underdense, in an otherwise homogeneous and expanding 
background universe. These spherical models provide a key 
reference for understanding and interpreting more complex 
situations. As a result of the spherical symmetry, the prob- 
lem is essentially one-dimensional, allowing a fully analytic 
treatment and solution, making the model easier to ana- 
lyze, interpret and understand. The spherical model for the 
evolution of isolated voids is discussed in some detail in Ap- 
pendix A. 
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Figure 3. Spherical model for the evolution of voids. Left: a pure (uncompensated) tophat void evolving up to the epoch of shell- 
crossing. Initial (linearly extrapolated) density deficit was Ai in q = —10.0, initial (comoving) radius Rio = 5.0h~ 1 Mpc. Timesteps: 
a = 0.05, 0.1, 0.2 and 0.3. Right: a void with an angular averaged SCDM profile (BBKS, cqn. 7.10). Initial density deficit and characteristic 
radius arc same as for the tophat void (left). The tendency of this void to evolve into a tophat configuration by the time of shell crossing 
is clear. Shell-crossing, and the formation of an obvious ridge, happens only if the initial profile is sufficiently steep. 



The most basic and universal properties of evolving 
spherical voids are: 

• Expansion: Voids expand, in contrast to overdense re- 
gions, which collapse. 

• Evacuation: As they expand, the density within them 
decreases continuously. (To first order, the density decrease 
is a consequence of the redistribution of mass over the ex- 
panding volume. Density decrease from mass lost to the sur- 
rounding overdensities is a second order, nonlinear effect oc- 
curing only near the edges.) 

• Spherical shape: Outward expansion makes voids evolve 
towards a spherical geometry. 

• Tophat density profile: The effective "repulsion" of the 
matter interior to the void decreases with distance from the 
center, so the matter distribution evolves into a (reverse) 
"tophat" . 

• "Super-Hubble" velocity field: Consistent with its (ul- 
timate) homogeneous interior density distribution, the (pe- 
culiar) velocity field in voids has a constant "Hubble-like" 
interior velocity divergence. Thus, voids evolve into genuine 
"Super-Hubble Bubbles". 

• Suppressed structure growth: Density inhomogeneities 
in the interior are suppressed and, as the object begins to 
resemble an underdense universe, structure formation within 
it gets frozen-in. 

• Boundary ridge: As matter from the interior accumu- 
lates near the boundary, a ridge develops around the void. 

• Shell- crossing : The transition from a quasi-linear to- 
wards a mature non-linear stage which occurs as inner shells 
pass across outer shells. 

Figure 3 illustrates these features. Both panels show the 
time evolution of the density deficit profile. Consider the 



panel on the left, which illustrates the development of an ini- 
tial (uncompensated) tophat depression (a "tophat" void). 
The initial (linear) density deficit of the tophat was set 
to Aii„,o = —10, and its (comoving) initial radius was 
Ri,o = 5/i _1 Mpc. The evolving density profile bears out 
the charactertistic tendency of voids to expand, with mass 
streaming out from the interior, and hence for the density 
to continuously decrease in value (and approach emptiness, 
8 — —1.0). Initially underdense regions are just expand- 
ing faster than the background and will never collapse (in 
an Q < 1 Universe). Notice that this model provides the 
most straightforward illustration of the formation of a ridge. 
Despite the absence of any such feature initially, the void 
clearly builds up a dense and compact bounding "wall" . 

For comparison with the tophat void configuration on 
the left, the panel on the right of Figure 3 depicts the evolu- 
tion of a void whose initial configuration is more representa- 
tive of cosmological circumstances. Here, the initial profile 
is the radial- averaged density profile for a trough in a Gaus- 
sian random field of Cold Dark Matter density fluctuations. 
The analytical expression for this profile was worked out by 
BBKS (1986) (eq. 7.10), and the one example we show here 
concerns the radial profile for a density dip with average 
steepness, i.e. x = — (V 2 /}/°"2 = — 1- The same qualita- 
tive aspects of void evolution can be recognized as in the 
case of a pure tophat void: the void expands, empties (to a 
near-empty configuration A = —1 at the centre), and also 
develops a ridge at its boundary. Notice that the void pro- 
file evolves into a configuration which increasingly resembles 
that of a "tophat" void. We will make use of this generic 
evolution in what follows. 

Looking from the inside out, one sees the interior shells 
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expanding outward more rapidly than the outer shells. With 
a minimum density near the void's centre, and density which 
increases gradually as one moves outward, the density deficit 
|A(r)| of the void decreases as a function of radius r. The 
outward directed peculiar acceleration is directly propor- 
tional to the integrated density deficit A(r, t) and therefore 
decreases with radius: inner shells are propelled outward at 
a higher rate so that the interior layers of the void move out- 
ward more rapidly. The inner matter starts to catch up with 
the outer shells, leading to a steepening of the density profile 
in the outer realms. Meanwhile, over a growing area of the 
void interior, the density distribution is rapidly flattening. 
This is a direct consequence of the outward expansion of the 
inner void layers: the "flat" part of the density distribution 
in the immediate vicinity of the dip gets "inflated" along 
with the void expansion. 

The features summarized above, which are seen in the 
idealized setting of initially smooth spherically symmetric 
voids, are also seen in more generic, less symmetric cos- 
mological circumstances, when substructure is also present. 
Figure 2 provides one illustration of the evolution of more 
realistic an complex underdensities. N-body simulation stud- 
ies of objects like this one have concluded that the tophat 
spherical model represent a remarkably succesfull desciption 
of reality (e.g. Dubinski et al. 1993; Van de Weygaert & Van 
Kampen 1993). The evolution towards a spherical top hat, 
whatever the initial configuration, is in stark contrast to 
how overdensities evolve. As a generic overdensity collapses, 
it contracts along a sequence of increasingly anisotropic con- 
figurations. Contraction leads to a "deflation" and accompa- 
nying steepening of density gradients, while the infall of sur- 
rounding structures marks a decreasing domain over which 
the neglect of substructure is realistic. 

In summary, it is apparent that the top hat spherical 
model not only provides a rather useful model for the evo- 
lution of isolated voids, but that it develops into an increas- 
ingly accurate representation of reality over an increasingly 
large fraction of the expanding volume of the void. 



3 EFFECT OF LARGER-SCALE STRUCTURES 

If we wish to use voids to understand the complex spatial 
patterns in the universe, we need a prescription for identi- 
fying the present-day cosmic voids, and for describing how 
voids interact with the large-scale structure which surrounds 
them. 



3.1 Importance of shell-crossing 

The generic property of ridge formation (e.g. Fig. 3) is sug- 
gestive, and Blumenthal et al. (1992) argued that the ob- 
served voids in the galaxy distribution should be identified 
with primordial underdensities that have only just reached 
shell- crossing. For a perfectly spherical void with a perfect 
tophat profile this happens exactly when the primordial den- 
sity depression out of which the void developed would have 
reached a linearly extrapolated underdensity <5 V . For such 
voids, 8 V is independent of mass scale: S v — —2.81 in an 
0,0 — 1 Universe. This threshold value will play an impor- 
tant role in our model of how the void hierarchy evolves. 



For instance, Dubinski et al. (1993) used this characteris- 
tic density to estimate that shell- crossing voids constitute 
a population of approximately volume- filling domains for a 
substantial range of cosmological structure formation sce- 
narios. By contrast, overdense primordial perturbations col- 
lapse and virialize — they shrink in comoving coordinates. 
The resulting picture is one in which the matter in the Uni- 
verse accumulates in ever smaller collapsing overdensities - 
in sheets, filaments and clusters - whose spatial arrangement 
is dictated by the growing underdense expanses. 



3.2 Void Sociology 

Two effects will seriously affect the number of small voids 
within a generic field of density perturbations. Both relate 
to the hierarchical embedding of a density depression within 
the larger scale environment. 

First, consider a small region which was less dense than 
the critical 5 V . It may be that this region, which we would 
like identify as a void today, was embedded in a significantly 
larger underdense region which was also less dense than the 
critical density. Therefore, we would also like to identify the 
larger region as a large void today. Since many small voids 
may coexist within one larger void, we must not count all 
of the smaller voids as distinct objects, lest we overestimate 
the number of small voids, and the total volume fraction 
in voids. We will call this the void-in-void problem. It is 
analogous to the well-known cloud-in- cloud problem associ- 
ated with the using the number density of initially overdense 
peaks to estimate the number of dense virialized clusters. 

A second effect is responsible for a radical dissimilarity 
between void and halo populations: If a small scale minimum 
is embedded in a sufficiently high large scale maximum, then 
the collapse of the larger surrounding region will eventually 
squeeze the underdense region it surrounds; the small-scale 
void will vanish when the region around it has collapsed 
completely. If the void within the contracting overdensity 
has been squeezed to vanishingly small size it should no 
longer be counted as a void. Figure 4 shows three examples 
of this process, each identified from a large (SCDM) N-body 
simulation. To account for the impact of voids disappearing 
when embedded in collapsing regions, we must also deal with 
the void-in-cloud problem. 

Virialized halos within voids are not likely to be torn 
apart as the void expands around them. Thus, the cloud- 
in-void phenomenon is irrelevant for dark halo formation. 
The asymmetry between the void-in-cloud and cloud-in-void 
processes effects a symmetry breaking between the emerging 
halo and void populations: although they evolve out of the 
same symmetric Gaussian initial conditions, we argue that 
over- and underdensities are expected to evolve naturally 
into agglomerations with rather different characteristics. 



4 EXCURSION FORMALISM 

In its simplest and most transparent formulation the excur- 
sion set formalism refers to the collapse of perfectly spherical 
overdensities, so this is the case which we will describe first. 
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Figure 4. Three examples (left to right) of the void-in- cloud process in action in numerical simulations of structure formation in a 
SCDM scenario (f2o = 1.0, ho = 0.5). Top to bottom panels show the evolution of the particle distribution in comoving coordinates from 
early to late times (resp. a = 0.3, 0.4, 0.5 and 1.0, the current epoch). The initially underdense regions are crushed by the collapse of the 
ovcrdense regions which surround them. The void in the first row of 4 panels shows a nearly spherical collapse sequence. The other two 
rows involve configurations involving more anisotropic surrounding matter distributions (and force fields). 



4.1 Excursion set model of clusters 

The jagged line in Fig. 5 represents the overdensity centred 
on a randomly chosen position in the initial Gaussian ran- 
dom field, as a function of the scale on which the overdensity 
was computed. The height of the walk S (S) is the linear 
theory overdensity relative to the density of the background 
universe. The spatial scale is parametrized by its variance S 
(defined in equation Bl). In hierarchical models, S decreases 
with increasing scale, so the largest spatial scales are on the 



left, and S(S) — > as S — > 0. Because the initial fluctua- 
tions are small, the mass contained within the smoothing 
filter is m oc [1 + (Di/ Do)5o] R 3 , where Di denotes linear 
theory growth factor at the initial time. Since D;/Do <JC 1, 
m oc R 3 : the mass is proportional to the initial comoving 
scale cubed. 

In the spherical collapse model, all regions with linear 
theory densities greater than <5 C can have formed bound viri- 
alized objects, and this critical overdensity is independent 
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Figure 5. Excursion Set Formalism, illustrated for the forma- 
tion of a halo. Random walk exhibited by the average overdcnsity 
8 centred on a randomly chosen position in a Gaussian random 
field, as a function of smoothing scale, parametrized by S m (large 
volume are on the left, small volumes on the right). Dashed hor- 
izontal line indicates the collapse barrier (5 C . The largest scale 
(smallest value of S) on which S(S) exceeds 8 C is an estimate of 
the mass of the halo which will form around that region. 



of mass scale. This constant value is shown as the dashed 
line in same height at all S m , where we have used the sub- 
script m to denote the fact that mass and initial scale are 
interchangeable. . 

The excursion set formalism supposes that no mass can 
escape from a region which collapses. If So — 8 C on scale 
R, then all the mass contained within R is included in the 
collapsed object, even if So < S c for all r < R. Thus, if the 
random walk height So exceeds the value <5 C after having trav- 
elled distance S(R) it represents a collapsed object of mass 
m oc R 3 . A walk may cross the barrier S c at many differ- 
ent values of S(R). Each crossing corresponds to a different 
smoothing scale and, because m oc _R 3 , contains a different 
amount of mass. However, of the various crossings of the 
barrier <5 C the first crossing, at the smallest value of S(R) 
for which So > 8 C , is special since it is this scale which is as- 
sociated with the most mass. The crossings at smaller scales 
correspond to condensations of a smaller mass, which have 
been incorporated in the larger encompassing mass concen- 
tration. 

In its simplest form, the excursion model for the dis- 
tribution of masses of virialized objects equates the distri- 
bution of distances S(R) which one-dimensional Brownian 
motion random walks, originating at the origin, travel be- 
fore they first cross a barrier of constant height S c , with the 
fraction of mass which is bound up in objects of mass m(R). 
The further a given walk travels before crossing the barrier, 
the smaller the mass of the object with which it is associated 
(Bond et al. 1991). 



4.2 Excursion set model of voids 

In our discussion above of the halo mass function, we consid- 
ered the cloud-in- cloud problem, and argued that the only 
cloud which should be counted was the largest possible one. 
To study voids in the excursion set approach one must first 
specify the boundary shape associated with the emergence of 
a void. This can be done if we know the critical underdensity 
5 V which defines a void, and in what follows we will use the 
epoch of shell-crossing, estimated using the spherical evolu- 
tion model, to specify 5 V . Thus, S v — —2.81, independent of 
smoothing scale (as was S c ). 

One might have thought that whereas clusters form 
from overdensities, voids form from underdensities, so the 
distribution of voids can be estimated analogously to how 
one estimates the distribution of clusters — one simply re- 
places the barrier S c with one at 8 V , and then studies the dis- 
tribution of first crossings of S v . Thus, if the random walk So 
first drops below the value 8 V after having travelled distance 
S(R) it represents a void of mass m cc R A and physical size 

n » i.7R. 

However, we have seen that we must be more careful; in 
addition to avoiding the double counting associated with the 
void-in-void process, we must also account for the void-in- 
cloud process. The strength of the excursion set formulation 
is that is shows clearly how to do this. Figure 6 illustrates 
the argument. There are four sets of panels. The left-most of 
each set shows the random walk associated with the initial 
particle distribution. The two other panels show how the 
same particles are distributed at two later times. The first set 
illustrates the cloud-in-cloud process. The mass which makes 
up the final object (far right) is given by finding that scale 
within which the linear theory variance has value S = 0.55. 
This mass came from the mergers of the smaller clumps, 
which themselves had formed at earlier times (centre panel). 
If we were to center the random walk path on one of these 
small clumps, it would cross the higher barrier S c /D(t) > S c 
at S > 0.55, the value of D(t) representing the linear theory 
growth factor at the earlier time t. 

The second series of panels shows the cloud-in- void pro- 
cess. Here, a low mass clump (S > 0.85) virializes at some 
early time. This clump is embedded in a region which is 
destined to become a void. The larger void region around it 
actually becomes a bona-fide void only at the present time, 
at which time it contains significantly more mass (S = 0.4) 
than is contained in the low mass clump at its centre. Notice 
that the cloud within the void was not destroyed by the for- 
mation of the void; indeed, its mass increased slightly from 
S > 0.85 to S ~ 0.85. Such a random walk is a bona-fide 
representative of S ~ 0.85 halos; for estimating halo abun- 
dances, the presence of a barrier at S v is irrelevant. On the 
other hand, walks such as this one allow us to make some 
important inferences about the properties of void-galaxies, 
which we will discuss shortly. 

The third series of panels shows the formation of a large 
void by the mergers of smaller voids: the void-in-void pro- 
cess. The associated random walk looks very much the in- 
verse of that for the cloud-in-cloud process associated with 
halo mergers. The associated random walk shows that the 
void contains more mass at the present time (S ~ 0.4) than 
it did in the past (5 > 0.4); it is a bona-fide representative of 
voids of mass S ~ 0.4. A random walk path centered on one 
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Figure 6. Four mode (extended) excursion set formalism. Each row illustrates one of the four basic modes of hierarchical clustering: 
the cloud- in- cloud process, cloud-in-void process, void-in-void process and void-in-cloud process (from top to bottom). Each mode is 
illustrated using three frames. Leftmost panels show 'random walks': the local density perturbation <5o(x) as a function of (mass) resolution 
scale Sm (cf. Fig. 5) at an early time in an N-body simulation of cosmic structure formation. In each graph, the dashed horizontal lines 
indicate the collapse barrier <5 C and the shell-crossing void barrier <5 V . The two frames on the right show how the associated particle 
distribution evolves. Whereas halos within voids may be observable (second row depicts a halo within a larger void), voids within collapsed 
halos are not (last row depicts a small void which will be squeezed to small size as the surrounding halo collapses). It is this fact which 
makes the calculation of void sizes qualitatively different from that usually used to estimate the mass function of collapsed halos. 
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of these mass elements which make up the filaments within 
the large void would resemble the cloud-in- void walk shown 
in the second series of panels. [Note that the height of the 
barrier associated with voids which are identified at cosmic 
epoch t scales similarly to the barrier height associated with 
halo formation: 8 v (t) = S v /D(t).] 

Finally, the fourth series of panels illustrates the void- 
in-cloud process. The particle distribution shows a relatively 
large void at the early time being squeezed to a much smaller 
size as the ring of objects around it collapses. A simple in- 
version of the cloud-in-void argument would have tempted 
one to count the void as a relatively large object containing 
mass S ~ 1. That this is incorrect can be seen from the fact 
that, if we were counting halos, we would have counted this 
as a cloud containing significantly more mass (S ~ 0.3), and 
it does not make sense for a massive virialized halo to host 
a large void inside. 

Thus, the excursion set model for voids which we will 
develop below is as follows: If a walk first crosses S c and 
then crosses 8 V on a smaller scale, then the smaller void is 
contained within a larger collapsed region. Since the larger 
region has collapsed, the smaller void within it no longer 
exists, so it should not be counted. The only bona-fide voids 
are those associated with walks which cross S v without first 
crossing <5 C . The problem of estimating the fraction of mass 
in voids reduces to estimating the fraction of random walks 
which first crossed 8 V at S, and which did not cross 8 C at any 
S' < S. Thus, a description of the void hierarchy requires 
solution of a two-barrier problem. 

Clearly, the model predictions will depend on 8 C and 
<5 V . If we use the spherical tophat model summarized in Ap- 
pendix A to set these values, then it seems reasonable to set 
<5 V = —2.81. But how we account for the void-in-cloud prob- 
lem is somewhat more subtle. Suppose we choose <5 C = 1.686, 
the value associated with complete collapse. In effect, this 
allows a void to have the maximum possible size it can have, 
given its underdensity, unless it is within a fully collapsed 
halo, in which case it has zero size. Presumably, if it is within 
a collapsing region which has not yet collapsed completely 
(as in the bottom-right panel of Figure 6), then its size is in- 
termediate between the size one would have estimated from 
the isolated spherical evolution model, and zero. Thus, only 
excluding voids in regions which have collapsed completely 
almost certainly overestimates the typical void size (further- 
more, we are ignoring the thickness of the ridge around each 
void). Another natural choice is <5 ta = 1.06; this ignores 
all voids that are within regions which are beginning to 
turnaround, even though they may still have non-neglibigle 
sizes, and so underestimates the abundance of large voids. 
Accounting more carefully for the effect of the void-in-cloud 
problem is the subject of ongoing work. 

In summary, what distinguishes voids from collapsed 
objects is the following: Whereas it may be possible to have 
a cluster within a void it does not make physical sense to 
have a void within a cluster. The excursion set formulation 
allows one to account for this. 



5 UNIVERSAL VOID SIZE DISTRIBUTION 

Let F(S, 8 V ,8 C ) denote the fraction of walks which first cross 
8 V at S, and which do not cross 8 C until after they have 



crossed 8 V (i.e., if they cross 8 C , they do so at s > S). 
Then J-(S, S v , S c ) is the distribution of first crossings of the 
type associated with voids. Appendix C shows that this first 
crossing distribution is given by 



SF(S,8 V ,8 C ) 



j 2 TY 2 V' 2 sin(jnV) 



x exp 



JIT 

" 2 8 2 /S 



where 



V = 



l«5v| 



(*c + |*v|) 



(1) 



(2) 



In an Einstein de-Sitter universe, S c , S v and a(m) all have 
the same time dependence, so equation (1) evolves self- 
similarly. In more general world-models the time depen- 
dences only slightly different, so the approximation of self- 
similar evolution should be quite accurate. 

The quantity V is the "void- and- cloud parameter"; it 
parameterizes the impact of halo evolution on the evolving 
population of voids. To see why, notice that the likelihood of 
smaller voids being crushed through the void-in-cloud pro- 
cess decreases as the relative value of the collapse barrier S c 
with respect to the void barrier 8 V becomes larger. 

This is also consistent with the fact that 



/ 



F{S,6 v ,6 e ) = l-V = 



8 C + \8 V 



(3) 



(e.g., equation C8) represents the mass fraction in voids. 
Thus, if V is small, voids account for nearly all the mass. 
On the other hand, for any noticeable impact of the void-in- 
cloud process the mass fraction in voids, 1 —T>, will be less 
than unity. The more important the void-in-cloud process 
is, the smaller the mass fraction in voids will be, as more 
voids are squeezed to vanishingly small size. 

Relation (7) suggests that the volume fraction in voids 
is 1.7 3 (1 - V). For S v = -2.81 and S c = 1.686, this ratio is 
larger than unity, indicating that the voids fill the universe. 
(The volume fraction in voids is also larger than unity if 
we set 8 C — 1.06 instead.) Thus, we have a model in which 
about one third of the mass of the universe is associated 
with voids which occupy most of the volume. The remaining 
seventy percent of the mass is in between the voids, and 
occupies negligible volume. 

Although the sum in equation (1) converges reasonably 
rapidly, it is not so easy to see what shape it implies. We 
have found that equation (1) is quite well approximated by 



"/M - exp {- 2 ) exp {-— l^- 2 -)' 



where we have set 



8 2 V /S. 



8 v /a (m) . 



(4) 



(5) 



and vf(y)Av/v — SJ-(S)dS/S. (This expression is accu- 
rate for values of 5 C / 1 <5 V | > 1/4 or so.) Expression (4) shows 
clearly that f(y) cuts-off sharply at both small and large 
values of v. In other words, the distribution of void masses 
is reasonably well peaked about v « 1, corresponding to a 
characteristic mass of order ao(m) £s <5 V |. 

When 8 C >• \8 V \, then V — > 0, and the second exponen- 
tial tends to unity. In this limit, the two-barrier distribution 
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Figure 7. Scaled distribution of void masses/sizes: voids which 
enclose large masses have large values of v. Curves show equa- 
tion (1) with <5 V = —2.81 and two choices of <5 C as labeled. These 
choices are motivated by the spherical collapse model, and re- 
sult in a distribution which is well peaked about a characteristic 
value. Ignoring the void-in-cloud process altogether is equivalent 
to setting S c — > oo. Although decreasing <5 C /|<5 V | decreases the 
abundance of low mass halos, the abundances of voids which en- 
close the most mass are not sensitive to the value of <5 C . 



reduces to that associated with a single barrier at S v . This 
shows explicitly that when the void-in-cloud process is unim- 
portant (T> — > 0), then the abundance of voids is given by 
accounting correctly for the void-in-void process. 

Figure 7 illustrates the resulting void size distributions. 
Notice that the mass fraction in small voids depends strongly 
on S c (the divergence at low v associated with the void-in- 
void solution is removed as S c increases), whereas the mass 
fraction enclosed by the largest voids depends only on S v . 
This is primarily a consequence of the fact that large un- 
derdensities embedded in a larger region of average density 
are rare, so such regions embedded in large overdensities are 
rarer still. Since there are essentially no large-scale under- 
densities embedded in larger scale overdensities, on scales 
where a <C (8 C + \ S V \), the value of <5 C is irrelevant. Thus, the 
distribution of large voids is almost exclusively determined 
by S v . We will return to this shortly. 

The number density n(m) of voids which contain mass 
m is obtained by inserting expression (1) in the relation 



m 2 n v (m) 



SF(S,5 V ,5 C ) 



d In S 
d km 



(6) 



To illustrate what our two-barrier model implies for void 
sizes, we must convert the expression above for the fraction 
of mass in voids to a void-size distribution. The simplest ap- 
proximation, motivated by the spherical tophat void model, 
sets the comoving volume v of the void equal to 



Since all the time dependence enters via v = S 2 (z)/a 2 (m), 
the distribution of void sizes evolves self-similarly. Simple 
changes of variables relate the void volume or mass functions 
to the barrier crossing distribution: f(v)dv = f(m)dm = 

As a specific illustration of what our model implies, 
Fig. 8 shows the distribution of void sizes in a model where 
the initial power-spectrum was P(k) oc k n with n = —1.5, 
normalized so that the rms fluctuations in a tophat sphere 
of radius unity was as — 0.9 at z — 0. The top left panel 
shows the mass fraction in voids of radius r, and the bot- 
tom left panel shows the number density of such voids. The 
three curves in each panel show equation (4) with S c = 1.06, 
1.686 and oo, and we have set S v — —2.81 in all cases. Notice 
how the abundance of small voids decreases dramatically as 
the ratio <J c /|<5v| decreases. By contrast, the abundance of 
large-scale voids is largely insensitive to this ratio (also see 
Fig. 7). 

We can make a rough estimate of the scale of the peak 
by computing that v at which equation (4) is maximized. 
This requires solution of a cubic, and gives i/ max decreasing 
as T> decreases. For the range of 0.6 < T> < 0.75 of interest, 
it is usually close to unity: i/ max w 1. To estimate the typical 
void size we will therefore simply use the approximate value 
of a ~ |<5 V |. 

For a power spectrum approximated by a power-law of 
slope n, the initial comoving size Ri of a region which is 
identifed as a void is, 



8 \(n+3)/2 



(I) 



Ri = 8 



as_ 

\5 V \ 



2/(n+3) 



(8) 



with ag denoting the rms fluctuation on scales of 8/i -1 Mpc 
(currently favoured ACDM models have as ~ 0.9). This 
means that the final size r v of the void is 



1.7 x 



8 



3 2/(3+n) 



as 2.7 
0.9 |5 V 



2/(3+n) 



(9) 



v = (m/p) x 1.7 3 



(7) 



A reasonable approximation to CDM spectra on Megaparsec 
scales is obtained by setting n = —1.5. In this case, the typ- 
ical void radius is ~ 3/i _1 Mpc. Since the correlation length 
is of order 8/i _1 Mpc, this makes the typical void diameter 
of order the correlation length. 

The top right panel shows the cumulative distribu- 
tion of the volume fraction for the three choices of S c . In 
all three cases, voids with radii greater than 5/i -1 account 
for about sixty percent of the volume. This suggests that, 
for sufficiently large voids, the details of the void-in-cloud 
process are not important. It is easy to see why: a typi- 
cal cluster forms from a region which had comoving radius 
Ri ~ 8(cr 8 /(5c) 2/(3+ ™ ) ~ 3.5/i _1 Mpc. Since few collapsing 
regions are larger than this, voids which are initially larger 
than this are extremely unlikely to have been squeezed out 
of existence. 

Finally we turn to an estimate of how the volume frac- 
tion in voids evolves in this model. Since as(z) = <7g/(l + z), 
the typical comoving size of voids is expected to be smaller 
at higher redshifts, by a factor of (1 + z) _2/(3+n) . The bot- 
tom panel shows the cumulative distribution at redshifts 
zero, one-half and unity (solid, dotted and dashed curves) 
where we have approximated S c (z) = 1.686(1 + z) and 
S v (z) = -2.81(1 + ^). 
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Figure 8. Distribution of void radii predicted by equation (4), in an Einstein de-Sitter model with P{k) oc k~ l s ; normalized to erg 
at 2 = 0. Top left panel shows the mass fraction in voids of radius r. Bottom left panel shows the number density of voids of radius 
r. Note that the void-size distribution is well peaked about a characteristic size provided one accounts for the void-in-cloud process. 
Top right panel shows the cumulative distribution of the void volume fraction. Dashed and solid curves in the top panels and bottom 
left panel show the two natural choices for the importance of the void-in-cloud process discussed in the text: <5 C = 1.06 and 1.686, with 
5 V = —2.81. Dotted curve shows the result of ignoring the void-in-cloud process entirely. Clearly, the number of small voids decreases 
as the ratio of <5 C /|<5 V | decreases. Bottom right panel shows the evolution of the cumulative void volume fraction distribution. The three 
curves in this panel are for <5 C = 1.686(1 + 2), where 2 = (solid), 0.5 (dotted) and 1 (dashed). 



: 0.9 



5.1 Alternative Models 

To better appreciate the ramifications of the two-barrier ex- 
cursion set model, it is instructive to explore alternative 
descriptions. This section discusses two models which follow 
from associating present-day voids with sufficiently under- 
dense troughs in the initial fluctuation field. 



5.1,1 The Basic Troughs Model 

The most straightforward model of the void distribution is 
to suppose that voids are associated with minima in the ini- 



tial density field. The simplest approximation to the number 
density of voids comes from smoothing the initial density 
fluctuation field with a filter of scale R, and then counting 
the number of minima of depth S v in the smoothed field. 
If one assumes that all the initial minima survive to the 
present time, then the number density of minima gives the 
number density of voids. BBKS (1986) show that the density 
of minima of depth 



(10) 



<T (m) 



in a Gaussian random field is 



© 0000 RAS, MNRAS 000, 000-000 



Void Hierarchy 13 



n v (u) Av 



u~exp(-v/2) G(7,7i/ 1/2 ) Av 



2ir (27r) 3 / 2 i?3 



(11) 



where the spectral parameters J?* and 7 depend on the 
shape of the power spectrum of the initial density fluctuation 
field, whose definition is given in Appendix B, along with 
that of the integral expression for the function G(7, 7^ 1,/2 ). 
(Strictly speaking, BBKS considered density maxima rather 
than minima. However Gaussian fluctuations are symmetric 
around the mean, so the density of peaks and troughs of the 
same absolute height is the same.) 

Notice that, in this model, the abundance of density 
minima in the primordial Universe depends on the depth of 
the minimum. If we define 

f(v) Av = (m/p) n v (v) Av 

and use the fact that the mass under a Gaussian filter is 

m = p(2iv) 3/2 R 3 , (12) 

then we have a quantity which one might interpret as the 
fraction of mass which is in minima of depth v. Unfortu- 
nately, for a comparison with the distribution of void sizes, 
this is a rather awkward quantity, since, in this picture, all 
voids contain the same mass m whatever their height v (be- 
cause the smoothing radius R is the same for all the voids). 

However, intuitively one would expect that deeper pri- 
mordial minima should be identified with voids containing 
more mass, something which the above expression does not 
accomplish self-consistently. The model discussed in the next 
subsection attempts to account for the correlation between 
void mass and depth. 
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Figure 9. Scaled distribution of void masses/sizes (equation 1): 
voids which contain large masses have large values of v. Ignoring 
the void-in- cloud process entirely yields the upper most curve. 
The spherical evolution model suggests 5c ~ I'M/2; in this case, 
the void distribution is reasonably well peaked about a character- 
istic value. Lower and upper dotted lines which extend to small 
values of v show predictions derived from equations (11) and (13) 
of the peaks model. These predictions depend on the shape of the 
initial power spectrum: curves assume P(k) oc fc . 



5.1.2 An Adaptive Troughs Model 

If, instead, we smooth the initial density field with a range 
of filter sizes R, and identify voids with minima of depth 
<S v /<7n(m), then, because <to decreases as J? oc m ly/3 in- 
creases, we have a model in which voids which contain 
more mass are associated with deeper minima. Appel & 
Jones (1990) show how the changing smoothing scale mod- 
ifies equation (11). The abundance of voids one obtains by 
replacing the BBKS (1986) formula (our equation 11) with 
the one given by Appel & Jones (1990) is 



~eacp(-i//2) #(7,7i/ /2 ) R 2 a\(R) Am/m 



(13) 



2tt 3(R*/R) 3 71/Va a 2(#) Av j v 

in which we have set m/p = (2ir) 3 ' 2 R 3 , the relation between 
mass and filter radius for a Gaussian smoothing filter, 7 is 
defined in equation (B7), and 

-y) 2 /2(l- 7 2 )] 



H{j,y) = / Ax x f{x) 



exp[ 



\/2tt(1-7 2 ) 



where f(x) is given in equation (B10), At large v (i.e., for 
deep minima), H ~ jv 1 ^ 2 G, where G is defined in Ap- 
pendix Bl, so this expression is the same as equation (11). 
The two expressions differ significantly at smaller v. If the 
initial spectrum of density fluctuations was a power law, 
P(k) oc k n , then equation (13) for the void mass function 
becomes 



"/(") = exp (--J —^375— 



5 + n\ 3 / 2 



(14) 



where we have used the fact that, for a Gaussian filter, 
7 2 = (n + 3)/(n + 5), and (R/R,) 2 = (n + 5)/6. Comparison 
with the excursion set approximation (equation 4) shows 
that both estimates contain the term y /, v/2n exp(— v/2), 
responsible for the exponential cut-off at large sizes. How- 
ever, the additional correction factors differ substantially. 
For instance, in contrast to the excursion set formula, the 
correction factor in this primordial troughs model explicitly 
depends on the shape of the initial power spectrum. 

Although the distribution of void sizes associated with 
equation (13) cuts off exponentially at large sizes, as does 
the excursion set formula, it diverges at small sizes: 



n(m) oc m 



(15) 



Since this peaks model ignores both the void-in-void and 
the void-in-cloud processes, the divergence towards small 
void sizes is likely to be a significant overestimate. However, 
the large scale cut-off is likely to be accurate, probably even 
more-so than the excursion-set approximation (see below). 

For comparison, the lower and upper dotted curves in 
Fig. 9 show the two predictions associated with these pri- 
mordial troughs models: equations (11) and (13). These 
predictions depend on the shape of the initial power spec- 
trum, and for the curves in Fig. 9 we have assumed P(k) oc 
k~ 15 . The contrast between the small scale divergence of 
the peak/troughs formulae, with the small scale cut-off 
for the excursion set distributions, is obvious. Notice that 
the peaks/troughs models predict systematically more very 
large voids than does the excursion set model. The reason 
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for this is closely related to the fact that the excursion set 
model does not include a factor like that in equation (8) 
of Appel & Jones (1990). For this reason, at large v, the 
peaks-based model is likely to be more accurate. 



5.2 Void distribution and spatial patterns 

Our extension of the excursion set formalism provides a use- 
ful framework within which to construct an understanding 
of the dichotomy between the overdense and underdense re- 
gions of space in any hierarchical structure formation sce- 
nario. 

Because voids occupy most of the volume, the peaked 
void distribution predicted by our excursion analysis has 
strong implications for the expected spatial patterns in the 
cosmic matter distribution. Since the sizes of most voids 
will be similar to the characteristic void size, our findings 
suggest that the cosmic matter distribution will resemble a 
foamlike packing of spherical voids of similar size and ex- 
cess expansion rate. The dynamical origin of such a mat- 
ter distribution had been recognized by various authors, in 
particular within the context of analyses based upon an ex- 
trapolation of the Zel'dovich approximation (see Shandarin 
& Zel'dovich 1989) and its extension, the adhesion approx- 
imation (Kofman, Pogosyan & Shandarin 1990). The role 
of voids in the latter had indeed been recognized by Sahni, 
Sathyaprakash & Shandarin (1994). These studies spurred 
the concepts of a cosmic web or cosmic skeleton (see e.g. 
Van de Weygaert 1991, Bond, Kofman & Pogosyan 1996, 
Novikov, Colombi & Dore 2003). Such patterns are naturally 
expected for cosmological scenarios with a lowpass power 
spectrum, characterized by a sharp spectral cutoff, as they 
would imply the imprint of an intrinsically dominating spa- 
tial scale. The four mode excursion formalism demonstrates 
and explains why the presence of such patterns is the nat- 
ural outcome for a considerably wider range of Gaussian 
structure formation models. 

As an interesting thought experiment, suppose we ex- 
trapolate our findings to an ultimate and asymptotic ex- 
treme: What if we approximate the "peaked" void distribu- 
tion by a "spiked" distribution centered on the characteristic 
void size? In such a scenario, the cosmic matter distribution 
would be organized by a population of equally sized, spher- 
ical voids, all expanding at the same rate, akin to the sce- 
nario suggested by Icke (1984). In this idealization, the walls 
and filaments would be found precisely at the midplanes be- 
tween expanding voids, and the resulting skeleton of the 
matter distribution would be precisely that of a Voronoi 
Tessellation (Voronoi 1908; Okabe et al. 2000 and references 
therein). Our results appear to offer an explanation for the 
fact that heuristic models, based upon the use of tessella- 
tions as spatial templates for the galaxy distribution, can 
succesfully reproduce a variety of galaxy clustering proper- 
ties (Van de Weygaert & Icke 1989; Van de Weygaert 1991; 
Goldwirth et al. 1995). 

6 THE VOID HIERARCHY 

The void distribution function derived in the previous sec- 
tion allows us to study in some detail the processes involved 



in the formation and development of void-dominated pat- 
terns in the cosmic matter distribution. We have already 
discussed such gross features as the void-filling factor, and 
the mass fraction in voids. But the excursion set analysis 
paves the way to a detailed assessment of the temporal de- 
pendence of a particular void, on its "ancestral" heritage as 
well as its spatial dependence on environmental factors. The 
following subsections touch upon a few of these elements of 
void evolution. 

6.1 Void mass and volume fractions 

We have already argued that, in an Einstein de-Sitter uni- 
verse, the mass fraction in voids does not evolve: approxi- 
mately one-third of the mass is in voids (equation 3), and 
that these voids fill space. This conclusion does not depend 
strongly on cosmological model. Because the collapse barrier 
5 v (a) decreases with time, the typical comoving void radius 
is larger at late times. Therefore, the mass contained within 
a typical void is larger at late times. On the other hand, 
the total mass fraction does not evolve, from which we infer 
that the small mass voids present at early times must merge 
with each other to make the more massive voids which are 
present at later times. 

6.2 Void ancestry 

The mass contained in a void at the present time was previ- 
ously partitioned up among many smaller voids, each sepa- 
rated by their own walls. This distribution can be estimated 
similarly to how Bond et al. (1991) and Lacey & Cole (1993) 
estimate the growth of clusters. 

Consider a void Vo which contains mass M at a time ao 
when the critical densities for spherical collapse turnaround 
and void shell- crossing are <5 c o and <5 v o, respectively. At an 
earlier epoch ai, the critical densities were 5 c i > <5 c o and 
\Svi\ > \S v o\- The fraction of M which was previously in 
voids that contained mass m at the earlier time a\ is given 
by inserting 

S c ^S cl -6 v0 ; 6 v ^6 vl -5 v0 ; S S(m) - S(M) (16) 

in equation (1). Integrating this over all possible ancestral 
voids (i.e., integrate over all < m < M), yields the mass 
fraction of M which was in voids at the earlier epoch also: 

/void(M) = 1+1 "J J— = J 7 • (!7) 

V del — Ovi / del — o v i 

Note how similar this expression is to the universal mass 
fraction in voids given by equation (3). Note in particular 
that this fraction is less than unity. This reflects the fact 
that, at earlier times, some of the mass currently affiliated 
with the void Vo was not part of the ancestral voids. In- 
stead, this fraction of its matter content resided in the walls 
(and filaments) which partitioned Vo into its many smaller 
constituent voids. In an Einstein-de-Sitter universe 

5 c i — Sco = 5 c0 (zi — zo), and 

Svi - Svo = 5 v o (zi - zo) , (18) 

so that the mass fraction of void matter which was in voids 
at the earlier time also is 

/void(M) = 1 - Vo^— ^r, (19) 
(l + 2l) 
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where T>o = <S v o |/(<5 c() — 5 v o) is the void- and- cloud parame- 
ter at the current epoch. Thus, at z\ « 20, this fraction is 
close to unity, whereas for large lookback times 21 3> 20 it 
tends to 1 — Do which is equal to the global void mass frac- 
tion (equation 3). In other words, the large voids emerging 
nowadays are to be traced back to an approximately average 
cosmic volume at early times. 

The transformations above allow one to write down the 
excursion set predictions for the rate at which smaller voids 
merge to make bigger ones. The calculation is analogous to 
the one used when estimating the merger rates of collapsed 
halos, and we will leave it for future work. In other words, 
one may reconstruct the ancestry of voids, the void merger 
tree, although this exercise will be complicated by the high 
rate of premature void mortality. 



quantify a fact which is intuitively obvious: that dense re- 
gions have a smaller fraction of their mass in voids. 
Furthermore, the typical void size scales as 

S(m) « S(M) + |<5 v -<5o(<5)| , (23) 

where the void size R(m) decreases as S(m) increases. Since 
|<5 V — So(8)\ increases as 8 increases, the typical void size is 
larger in regions of lower density. Moreover, the sharpness 
of the peak in the void size distribution depends on <5 C /|<5 V | 
(c.f. Fig. 9: the void size distribution becomes more more 
sharply peaked as void-in-cloud demolition becomes more 
important. The transformations in equation (21) mean that, 
in dense regions (8 > 0), where voids are more likely to 
be demolished by collapsing clouds, the distribution of void 
sizes is expected to be narrower. 



6.3 Environmental Dependence 



6.4 Spatial Clustering 



Suppose we evaluate the density field smoothed on a grid 
with cells of size R. The smoothed density will fluctuate 
from cell to cell. In the excursion set approach, we find that 
voids in denser cells 1) are smaller, 2) have a narrower size 
distribution, and 3) account for a smaller fraction of to- 
tal mass in the cell they inhabit. This subsection quantifies 
these trends of "void bias" . 

Consider a cell of size V within which the density is 
p(l + 8); i.e., this cell contains mass M = pV(l + 8). In the 
spherical evolution model, the initial and final densities are 
related: 



So(S) = 



Ssc 



1.68647 



1.68647 



1.35 



(l + 5) 2 /3 



1.12431 



+ 



0.78785 



(1+<5)V2 (1 + 5)0.58661 



(20) 



(e.g. Mo & White 1996). Note that 8 has the same sign as 8; 
initially dense regions become denser, whereas the comoving 
density in underdense regions decreases with time. 

In the context of the void model studied here, voids 
which are in cells of volume V within which the over- 
density is 8 are described by random walks which do not 
start from the origin [S = 0, 80 = 0], but from the posi- 
tion [S(M),5o(S)]. Therefore the fraction of the total mass 
M = pV(l + 8) which is in voids of mass m is given by 
setting 



8 c -8 (8), 



Sv — > S v — So(8) 



S(m)-S pV(l + S) 



(21) 



in equation (1). Integrating the resulting distribution over 
< m < M yields the fraction of mass, in a region of volume 
V within which the density is 8, which is contained in voids: 



/void (8) 



5 (6) 



(22) 



This indicates that the mass fraction / vo id(<5) decreases as 
the density 8 of the cell increases. Conversely, as S — > —0.8, 
the density we associate with a void, then 8o(S) — * 8 V , and 
so / V oid(<5) — > 1 as expected. (In this extreme, the fitting 
formula (20) is slightly inaccurate, since it sets 8o(— 0.8) = 
—2.7, rather than —2.81.) Thus, our analysis allows one to 



The model developed here also allows us to build an ap- 
proximate model of the evolution of the dark matter cor- 
relation function following methods outlined in Neyman & 
Scott (1952) and Scherrer & Bertschinger (1992) (recently 
reviewed by Cooray & Sheth 2002). The calculation requires 
estimates of 1) the distribution void sizes, 2) the clustering 
of void centres on large scales, and 3) the density run within 
a void. The previous sections derived estimates for the first 
of these three quantities. 

The second one, the clustering of void centres, can be 
estimated as follows. Write the two-point correlation func- 
tion of voids which contain mass mi and m2 as 



^vv(r|mi,m 2 ) = b{m\)b{m,2) fdm(r), 



(24) 



where £ dm is the correlation function of the dark matter, and 
the bias factor b(m) depends on the mass or size of the voids. 
Following Cole & Kaiser (1989), Mo & White (1996) and 
Sheth & Tormen (1999), knowledge of the number density of 
objects is sufficient for estimating their spatial distribution, 
at least on large scales. Therefore, b(m) depends on which 
estimate of n v (m) we use. If we use equation (13) from the 
peaks model, then 



b(m) = 1 + 



v_ + hi 

I'M ' 



where 



(-l) k ( 7 ^ 1/2 ) fc d k H^,x) 



k\ (7,71/1/2) dx k 



(25) 



(26) 



1,1/2 



Using our approximation to the excursion set prediction 
(equation 4) instead gives 



\Sy\ 4u(S c + \S v \)- 

In both cases, the largest voids are more strongly clustered 
than those of average size. The higher order moments of 
the void distribution can be estimated similarly to how Mo, 
Jing & White (1997) estimate the higher order moments of 
clusters. 

If we suppose that all the mass is contained in the void 
walls, then we can approximate the density run around a 
void centre as a uniform density shell. Figures 3 and 5 in 
Dubinski et al. (1993) suggest this is a fair approximation. 
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Specifying the mass associated with the void as well as the 
shell thickness sets the density within the shell. Thus, we 
have all three ingredients required to model the power spec- 
trum (or correlation function) of the dark matter distribu- 
tion. 

There is one important aspect in which this void-based 
model for the correlation function differs from the usual 
halo-based model. Namely, in the halo model, halos arc 
treated as hard spheres which do not overlap; this leads to 
exclusion effects on small scales. Since the radius of a typical 
collapsed halo is smaller than a Megaparsec, the effects of 
exclusion are expected to be unimportant. In a void-based 
model, on the other hand, typical void radii are of order a few 
Megaparsecs; since voids do not overlap, exclusion effects are 
likely to matter on scales of order a few Megaparsecs. We 
leave a more extensive analysis of all this to future work. 



7 SUMMARY AND INTERPRETATION 

Initially underdense regions expand faster than the Hubble 
flow. If they are not embedded within overdense regions, 
such regions eventually form voids which are surrounded by 
dense void walls. These voids expand with respect to the 
background Universe, and during their expansion tend to 
become more and more spherical (Figure 2). The outward 
expansion is differential, so most initial void configurations 
tend to evolve to distinct "tophat" density profiles (Fig- 
ure 3). A description of the evolution of initially spherical 
tophat over- and underdense regions has been available for 
some time (Appendix A). Although the spherical evolution 
model allows one to study the evolution of single isolated 
objects, a more complete theory must also describe void 
evolution within the context of a generic random density 
fluctuation field. 

The evolving void hierarchy is determined by two pro- 
cesses: 

• The void-in-void process describes the evolution of a 
system of voids which are embedded in a larger scale under- 
density; in this case small voids from an early epoch merge 
with one another to form a larger void at a later epoch 
(Fig. 2). 

• The void-in- cloud process is associated with underdense 
regions embedded within a larger overdense region; in this 
case the smaller voids from an earlier epoch may be squeezed 
out of existence as the overdense region around them col- 
lapses (Fig. 4). 

In contrast, the evolution of overdensities is governed only by 
the cloud-in- cloud process; the cloud-in-void process is much 
less important, because clouds which condense in a large 
scale void are not torn apart as their parent void expands 
around them. 

This asymmetry between how the surrounding environ- 
ment affects halo and void formation can be incorporated 
into the excursion set approach by using one barrier to model 
halo formation and a second barrier to model void formation 
(Fig. 6). Only the first barrier matters for halo formation, 
but both barriers play a role in determining the expected 
abundance of voids. The resulting void size distribution is 
a function of two parameters (equation 1), which the model 
associates with the dynamics of expansion and collapse. The 



predicted distribution of voids is well-peaked about a charac- 
teristic size (Figs. 7 and 8) — in contrast, the distribution of 
halo masses is not. Comparison of the two-parameter family 
of void distribution curves (Figure 9) with the void size dis- 
tribution in numerical simulations of hierarchical clustering 
is the subject of work in progress (Colberg et al. 2004). 

Five major observations about the properties of the void 
population result from the two-barrier excursion set model: 

• The void-in-cloud mechanism (Fig. 4) is responsible for 
the demise of a sizeable population of small voids. As a re- 
sult, the void size distribution has a small-scale cut-off: the 
void population is "void" of small voids (Section 5), in a way 
which our excursion set analysis quantifies. 

• The population of large voids is insensitive to this effect 
(Fig. 7). Therefore, the abundance of voids which are larger 
than the typical initial comoving sizes of clusters should be 
well described by peaks theory or its extensions described 
in Section 5.1. 

• At any cosmic epoch there is a characteristic void size 
which increases with time: the larger voids present at late 
times formed from mergers of smaller voids which formed at 
earlier times (e.g Fig. 2 and Section 6). 

• At any given time the mass fraction in voids is approx- 
imately thirty percent of the mass in the Universe, and the 
voids approximately fill space (Section 6). 

• As the size of most voids will be similar to the charac- 
teristic void size, the cosmic matter distribution resembles a 
foamlike packing of spherical voids of approximately similar 
size and excess expansion rate. This may explain why sim- 
ple models based on the Voronoi tesselation exhibit many 
of the features so readily visible in N-body simulations of 
hierarchical clustering. 

7.1 Galaxies in Voids 

It is with some justification that most observational atten- 
tion is directed to regions where most of the matter in the 
Universe has accumulated. Almost by definition they are the 
sites of most observational studies, and the ones that arc 
most outstanding in appearance. Yet, for an understanding 
of the formation of the large coherent foamlike patterns per- 
vading the Universe, it may be well worth directing atten- 
tion to the complementary evolution of underdense regions. 
These are the progenitors of the observed voids, the vast re- 
gions in the large-scale cosmic galaxy distribution that are 
practically devoid of luminous matter. 

When extensive systematic redshift surveys began map- 
ping the spatial galaxy distribution, voids were amongst 
the most visually striking features. Since then, the role of 
voids as key ingredients of the cosmic galaxy distribution 
has been demonstrated repeatedly in extensive galaxy red- 
shift surveys (see Kauffmann & Fairall 1991; El-Ad, Pi- 
ran & da Costa 1996; El-Ad & Piran 1997; Hoyle & Vo- 
geley 2002; Plionis & Basilakos 2002; Rojas et al. 2003). 
A number of studies also indicated that observed voids ex- 
hibit distinct hierarchical features. Van de Weygaert (1991) 
suggested the existence of a "void hierarchy" when point- 
ing out that the galaxy distribution in the CfA/SRSS2 red- 
shift survey (Geller & Huchra 1989; da Costa et al. 1993) 
gave the impression of small-scale voids embedded in the less 
pronounced large-scale underdense region delimited by the 
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"Great Wall" . Even in the most canonical specimen amongst 
its peers, the Bootes void, traces of a faint structured inter- 
nal galaxy distribution were found (Szomoru et al. 1996). 

The dynamical impact of voids has proven to be cru- 
cial for understanding the cosmic flow patterns in the Local 
Universe. Measured peculiar galaxy velocities imply recon- 
structions of the local cosmic density field in which the re- 
pulsive actions of voids are important (e.g. Bertschinger et 
al. 1990; Strauss & Willick 1995; Dekel & Rees 1995). And 
more locally, the void's influence on cosmic flows was estab- 
lished when Bothun et al. (1992) studied galactic peculiar 
motions along a wall around the largest void in the CfA 
redshift sample (de Lapparent, Geller & Huchra 1986). 

In all these respects, voids in the galaxy distribution 
are similar to those in the dark matter distribution. How- 
ever, although voids in the galaxy distribution are mostly 
roundish in shape, they have typical sizes in the range of 
20k -1 - 50/i _1 Mpc (e.g. Hoyle & Vogeley 2002; Plionis & 
Basilakos 2002; Arbabi-Bidgoli & Midler 2002). These sizes 
are considerably in excess of the typical void diameters in 
our model of voids in the dark matter distribution, but note 
that the typical void size in the galaxy distribution depends 
on the galaxies which were used to define the void. The 
voids associated with rare luminous galaxies are larger in 
part because the number density of such galaxies is lower. 
As we describe below, our excursion set analysis provides a 
framework for modeling this dependence. 

In recent years, the possibility that void galaxies are a 
systematically different population has received considerable 
attention (see e.g. Szomoru et al. 1996; El- Ad & Piran 2000; 
Peebles 2001; Mathis & White 2002, Rojas et al. 2003; Ben- 
son et al. 2003). In the simplest models of biased galaxy 
formation (e.g. Little & Weinberg 1994) one would expect 
to find voids filled with galaxies of low luminosity, or galax- 
ies of some other uncommon nature (e.g. Hoffman, Silk & 
Wyse 1992). Indeed, even though various studies were ori- 
ented towards establishing the properties of voids in galaxies 
(e.g., Kauffmann & Fairall 1991; El-Ad & Piran 1997, 2000; 
Hoyle & Vogeley 2002; Arbabi-Bidgoli & Miiller 2002; Plio- 
nis & Basilakos 2002), and some focussed explicitly on the 
identity of galaxies inside voids (e.g. Szomoru et al. 1996; 
El- Ad & Piran 2000; Rojas et al. 2003), a clear picture of 
the relation between void galaxies and their surroundings 
is only just becoming available. This is in large part due to 
the fact the large scale surveys such as the SDSS (Abazajian 
et al. 2002) and 2dFGRS (Colless et al. 2003) now probe a 
sufficiently large cosmological volume that they contain a 
statistically significant number of large voids. 

Recently, Mathis & White (2002) and Benson et 
al. (2003) have identified and studied voids and void galaxies 
in semi-analytic galaxy formation models. In these models, 
the properties of galaxies are determined by the halos they 
inhabit. Therefore, if one can model the halo population as- 
sociated with voids, a model of the void galaxy population 
is within reach. The excursion set model developed here is 
phrased in the same language used in the simulations, so 
it represents the ideal framework within which to attempt 
such a model. 

In particular, consider the cloud-in-void process shown 
in the second series of panels in Fig. 6. Notice that the con- 
dition that the cloud exist in a void means that, on average, 
clouds in voids will be less massive than clouds in regions of 



average density (to represent a cloud, the walk must reach 
8 C , and on average, it will take more steps to travel to 5 C from 
8 V than from zero — more steps imply smaller masses). For 
similar reasons, the clouds associated with the more mas- 
sive halos should be more massive on average (this is also 
discussed more fully by Mo & White 1996; Sheth & Tor- 
men 2002 and Gottlober et al. 2003). Although we speak of 
the clouds as being within the voids, our discussion of how 
voids empty their mass into the ridge which surrounds them 
(c.f. Fig. 3) suggests it may be more appropriate to think of 
these clouds as being associated with the void walls. It seems 
natural to associate void galaxies with such clouds-in-voids. 
If low mass halos host lower mass galaxies, and less massive 
galaxies tend to be less luminous and bluer, then void galax- 
ies should be fainter and bluer than field or cluster galaxies; 
our model allows one to quantify this trend. Thus, the re- 
sults presented here allow a more elaborate model for voids 
in the galaxy distribution and the galaxy population in voids 
than that discussed recently by Friedmann & Piran (2001). 
Developing such a model is the subject of work in progress. 



ACKNOWLEDGEMENTS 

We thank Jorg Colberg, Bernard Jones and Paul Schechter 
for encouraging discussions and suggestions, and A. Babul 
for useful conversations on collapsing voids. This work was 
supported in part by the DOE and NASA grant NAG 5- 
10842 at Fermilab. 

O, what men dare do! What men may do! What men daily 
do, not knowing what they do! (Shakespeare 1598) 



REFERENCES 

Abazajian K., et al. (The SDSS collaboration) 2003, in press 

(www. sdss . org) 
Aikio J., Mahonen P., 1998, ApJ, 497, 534 
Appel L., Jones B. J. T., 1990, MNRAS, 245, 522 
Bardeen J.M., Bond J.R., Kaiser N., Szalay A., (BBKS), 1986, 

ApJ, 304, 15 

Benson A.J., Hoyle F., Torres F., Vogeley M.S., 2003, MNRAS 
340, 160 

Bertschinger E., 1983, ApJ, 268, 17 
Bertschinger E., 1985, ApJS, 58, 1 

Bertschinger E., Dekel A., Faber S.M., Dressier A., Burstein D., 

1990, ApJ, 364, 370 
Blumcnthal G.R., Da Costa L., Goldwirth D.S., Lecar M., Piran 

T., 1992, ApJ, 388, 234 
Bond J. R., Cole S., Efstathiou C, Kaiser N., 1991, ApJ, 379, 

440 

Bond J.R., Kofman L., Pogosyan D., 1996, Nature, 380, 603 
Bond J. R., Myers S., 1996, ApJS, ??? 

Bothun CD., Geller M.J., Kurtz M.J., Huchra J. P., Schild R.E., 

1992, ApJ, 395, 347 

Colless et al. (2dF consortium), 2003, astro-ph/0306581 

Cooray A., Sheth R. K., 2002, Phys. Rep., 372, 1 

da Costa L.N. 1993, in Cosmic Velocity Fields, Proc. 9th IAP 

Astrophysics Meeting, eds. F.R. Bouchet & M. Lachieze-Rey 

(Editions Fronticres), p. 475 
Dekel A., Rees M., 1994, ApJ, 422, LI 

De Lapparent V., Geller M.J., Huchra J. P., 1986, ApJ, 302, LI 
Dubinski J., Da Costa L.N., Goldwirth D.S., Lecar M., Piran T., 

1993, ApJ, 410, 458 

El-Ad H., Piran T., da Costa L.N., 1996, ApJ 462, L13 



© 0000 RAS, MNRAS 000, 000-000 



18 R. K. Sheth & R. van de Weygaert 



El-Ad H., Piran T., 1997, ApJ, 491, 421 
El-Ad H., Piran T., 2000, MNRAS, 313, 553 
Epstein R. A., 1983, MNRAS, 205, 207 

Einasto J., Einasto M., Gramann M., 1989, MNRAS, 238, 155 
Eiscnstcin D., Locb A., 1995, ApJ, 439, 520 
Fillmore J. A., Goldreich P., 1984, ApJ, 281, 9 
Friedmann Y., Piran T., 2001, ApJ, 548, 1 
Geller M, Huchra J., 1989, Science, 246, 897 
Goldberg D.M., Vogeley M.S., 2003, astro-ph/0307191 
Goldwirth D.S., Da Costa L.N., Van de Weygaert R., 1995, MN- 
RAS, 275, 1185 

Gottlober S., Lokas E., Klypin A., Hoffman Y., 2003, MNRAS, 

subm., astro-ph/0305393 
Gunn J. E., Gott J. R., 1972, ApJ, 176, 1 

Hausman M.A., Olson D.W., Roth B.D., 1983, ApJ, 270, 351 

Heavens A., Peacock J., 1988, MNRAS, 232, 339 

Hoffman G.L., Salpcter E.E., Wasserman I., 1983, ApJ, 268, 527. 

Hoffman Y., Silk J., Wyse R. F. G., 1992, ApJ, 388, L13. 

Hoffman Y., Shaham J., 1982, ApJ, 262, L23 

Hoyle F., Vogeley M., 2002, ApJ, 566, 641 

Icke V., 1972, Formation of Galaxies Inside Clusters, Ph.D. The- 
sis, Leiden Univ. 
Icke V., 1973, A&A, 27, 1 
Icke V., 1984, MNRAS, 206, 1 

Kauffmann G., Fairall A. P., 1991, MNRAS, 248, 313 
Kauffman G., White S.D.M., 1993, MNRAS, 261, 921 
Kofman L., Pogosyan D.Yu., Shandarin S.F., 1990, MNRAS, 242, 
200 

Lacey C, Cole S., 1993, MNRAS, 262, 627 

Lilje P.B., Lahav O., 1991, ApJ, 374, 29 

Lin C. C, Mestel L., Shu F. H., 1965, ApJ, 142, 1431 

Little B., Weinberg D.H., 1994, MNRAS, 267, 605 

Mathis H., White S. D. M., 2002, MNRAS, 337, 1993 

Mo H. J., White S. D. M., 1996, MNRAS, 282, 347 

Mo H. J., Jing Y., White S. D. M., 1997, MNRAS, 284, 189 

Monaco P.L., 1997, MNRAS, 287, 753 

Ncyman J. Scott E. L., 1952, ApJ, 116, 144 

Novikov D., Colombi S., Dore, 2003, MNRAS, submitted, astro- 
ph/0307003 

Arbabi-Bidgoli, Muller V., 2002, MNRAS, 332, 205 

Peebles P.J.E., 1980, The Large Scale Structure in The Universe. 

Princeton Univ. Press, Princeton, NJ. 
Peebles P.J.E, 1982, ApJ, 257, 428 
Peebles P.J.E., 2001, ApJ, 557, 495 
Plionis M., Basilakos S., 2002, MNRAS, 330, 399 
Press W., Schechter P., 1974, ApJ, 187, 425 

Okabe, A., Boots, B., Sugihara, K., & Nok Chiu, S. 2000, Spa- 
tial Tessellations, Concepts and Applications of Voronoi Dia- 
grams, 2nd edition, Chichester, John Wiley & Sons Ltd 

Regos E., Geller M., 1991, ApJ, 377, 14 

Rojas R.R., Vogeley M.S., Hoyle F., Brinkmann J., 2003, astro- 

ph/0307274 
Ryden B., 1995, ApJ, 452, 25 

Sahni V., Sathyaprakash B.S., Shandarin S.F., 1994, ApJ, 431, 
20 

Scherrer R.J., Bertschinger E., 1991, ApJ, 381, 349 
Schucker P., Bohringer H., Arzner K., Reiprich T.H., 2001, A&A, 
370, 715 

Shakespeare W., 1598, Much Ado About Nothing 

Shandarin S.F., Zcl'dovich Ya.B., 1989, Rev. Mod. Phys., 61, 185 

Shectman S. A., Landy S.D., Ocmlcr A., Tucker D.L., Lin H., 

Kirshner R. P., Schechter P.L., 1996, ApJ, 470, 172 
Sheth R. K., 1998, MNRAS, 300, 1057 
Sheth R. K., Tormcn G., 1999, MNRAS, 308, 119 
Sheth R. K., Mo H., Tormen G., 2001, MNRAS, 323, 1 
Sheth R. K., Tormen G., 2002, MNRAS, 329, 61 
Suto Y., Sato K., Sato H., 1984, Prog. Theor. Phys., 71, 938 



Szomoru A., van Gorkom J., Gregg M.D., Strauss M.A., 1996, 

AJ, 111, 2150 
Strauss M., Willick J., 1995, Phys. Rep. 261, 271 
Van de Weygaert R., Icke V., 1989, A&A, 213, 1 
Van de Weygaert R., 1991, Voids and the Large Scale Structure 

of the Universe, PhD thesis, Leiden Univ. 
Van de Weygaert R., Bertschinger E., 1996, MNRAS, 281, 84 
Van de Weygaert R., Van Kampcn E., 1993, MNRAS, 263, 481 
Van de Weygaert R. , 2002, in Proceedings 2nd Hellenic Cosmol- 
ogy Workshop, Astrophysics and Space Science Library, vol. 
276, eds. M. Plionis & S. Cotsakis, Kluwer, pp. 119-259 
Voronoi G., 1908, J. reine angew. Math., 134, 198 
White S. D. M., 1979, MNRAS, 189, 831 
White S.D.M., Silk J., 1979, ApJ, 231, 1 
Zchavi E., et al., 2002, ApJ, 571, 172 

APPENDIX A: THE SPHERICAL TOPHAT 
MODEL 

Al Background 

Analytically tractable idealizations help to understand var- 
ious aspects of void evolution. In this regard, the spherical 
model represents the key reference model against which we 
may assess the evolution of more complex configurations. 
Also, it provides the clearest explanation for the various 
void characteristics listed in the main test. And most sig- 
nificantly within the context of this work, it provides the 
fundament from which our formalism for hierarchical void 
evolution is developed. 

The structure of a spherical void or peak can be treated 
in terms of mass shells. In the "spherical model" concentric 
shells remain concentric and are assumed to be perfectly 
uniform, without any substructure. The shells are supposed 
never to cross until the final singularity, a condition whose 
validity is determined by the initial density profile. The re- 
sulting solution of the equation of motion for each shell may 
cover the full nonlinear evolution of the perturbation, as long 
as shell crossing does not occur. 

The treatment of the spherical model in a cosmologi- 
cal context has been fully worked out (Gunn & Gott 1973; 
Lilje & Lahav 1991). As long as the mass shells do not cross, 
they behave as mini-Friedmann universes whose equation of 
motion assumes exactly the same form as that of an equiva- 
lent FRW universe with a modified value of Q a . The details 
of the distribution of the mass interior to the shell are of 
no direct relevance to the evolution of each individual shell. 
Instead, the evolution depends on the total mass contained 
within the radius of the shell, and the global cosmological 
background density. 

Although quantitative details depend on the cosmologi- 
cal model, a study of the evolution of spherical perturbations 
in an Einstein-de Sitter Universe suffices to illustrate all the 
important physical features. 



A2 Definitions 

When a mass shell at some initial time i; starts expanding 
from a physical radius r; = a(ti)xi, its subsequent motion is 
characterized by the expansion factor lZ(t,ri) of the shell: 

r(t,n) = n{t,ri)n, (Al) 
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where r(t, n) = a(t)x(t,Xi) is the physical radius of the shell 
at time t and x(t,Xi) the corresponding comoving radius. 
The evolution of the shell is dictated by the cosmological 
density parameter 



8TvGp u (t) 
3H?, 



(A2) 



and the mean density contrast within the radius of the shell, 



AM) = 



[ 


p(y,t) 1 




Pu{t) 



y 2 dy 



r 3 l 



S(y,t) y 2 dy, 



(A3) 



To determine the evolution of R(t,ri), it is convenient to 
introduce the parameters A C i = A c (ii) and en where 



1 + A ci = Hi [1 + A(ti,n) 



- 1. 



(A4) 



Here, fii = fi (ti), Hi = H{U) and Vi is the physical velocity 
(i.e. the sum of the peculiar velocity and Hubble expansion 
velocity with respect to the void center) of the mass shell 
at t = U. The usual assumption of a growing mode per- 
turbation implies that the velocity perturbation i> pec ,; for a 
spherical perturbation, at the initial time ti, is 

*w = -^p/tnOAfc.ti), (A5) 

and hence, 

on = -|/(ni)A(ri,ti). (A6) 

In effect, A c ; is the density contrast of the shell with respect 
to a critical universe (51 = 1) at the cosmic time ti, while 
at is a measure of the corresponding peculiar velocity (or, 
rather, the kinetic energy) of the shell. The evolution of a 
spherical over- or underdensity is entirely and solely deter- 
mined by the initial (effective) over- or underdensity within 
the (initial) radius r; of the shell, A c j(r,,tj), and the corre- 
sponding velocity perturbation, v pe c,i- Hence, the values of 
A c i and ai determine whether a shell will stop expanding or 
not, i.e. whether it is closed, critical or open. The criterion 
for a closed shell is A ci > ai, for a critical shell, A ci — at, 
and A c i < a; for an open shell. 

Notice that these expressions assume that the initial 
density fluctuation was negligible, so that the initial mass 
m and initial comoving size Ti are related: m oc R 3 . 



A3 Shell Solutions 

The solution for the expansion factor 1Z(t, n) = TZ(Q r ) of an 
overdense cq. underdense shell is given by the parametized 
expressions 



(A7) 



\ rl + Ac \ (l-cos8 r ) A ci >ai, 

yL±ci ai J 



in which the development angle O r , which paramaterizes all 
physical quantities relating to the mass shell, is related to 
time t via 



t(O r ) 



1 1 + A e 

2 (aj — A C i 



1 1 + A c 

2 (A C i — at 



(sinhO r — 9 r ) 



(6 r — sin Or 



A ci < a it 



(A8) 



i-V 2 



A c i > ai . 



while for a critical shell the solution is given by the direct 
relation 



2/3 



Aci = «i 



(A9) 



Notice that the solutions for the evolution of overdense and 
underdense regions in essence are the same, and are inter- 
changeable by replacing 



(sinhe - 9) 
(coshG- 1) 



(9 - sin 9) 

(1-COS0). 



(A10) 



A4 Density Evolution 



If the initial density contrast of a shell is Aj(ri), its density 
contrast A(r,t) at any subsequent time t is given by 



1 + AM) = 



1 + Ai(n) a(tf 



TI 3 



(All) 



With A(r, t) being a relative quantity, comparing the density 
of the mass shell at radius r at time t with that of the 
global cosmic background, the value of A(r, t) is a function 
of the shell's development angle 6 r as well as that of the 
development angle of the Universe Q u , 

2 



cosh O u + 1 
2 



fi < 1, 

a > i. 



n = { (A12) 

O ^ 1 

cos O u + 1 

The shell's density contrast may then be obtained from 

i + AM) = /(e r )//(e„), (Ai3) 

where /(G) is the cosmic "density" function: 



/(e) 



( (sinhQ-Q) 2 
(coshG - l) 3 

2/9 

(Q-sinQ) 2 

k (i - cose) 3 



open , 
critical , 

closed , 



(A14) 



This expression is equally valid for the shell (in which case 
"open" means A C i < ai) and the global background Uni- 
verse (where "open" means f2 < 1). 



A5 Shell Velocities 

The velocity of expansion or contraction of a spherical shell 
is given by computing d7?./dt, so it can be written in terms 
of 6 r and Q u . In particular, the shell's peculiar velocity with 
respect to the global Hubble velocity, 



v pec (r,t) = v(r,t) - H u (t)r(t) , 
may be inferred from the expression 

twM) = H u (t)r(t) j^g - 1 



(A15) 



(A16) 
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where H u (t)r = (a/a)r and the cosmic "velocity" function 



-D(Zcoll) Ahtj.O = 8 C 



(A23) 



3(6) = { 



' sinh 9 (sinh 6- 6) 
(cosh 6- l) 2 

2 

sin 6(6 -sin 6) 
k (l-cos6) 2 



open , 
critical , 
closed . 



(A17) 



Thus, we may define a Hubble parameter H s for each indi- 
vidual shell, 



(A18) 



A6 Overdensities and collapse when $7 = 1 

The previous sections provided explicit expressions for the 
evolution of a spherical perturbation in FRW backgrounds 
with no cosmological constant. To better illustrate our ar- 
gument, we will now specialize to the case of an Einstein 
de-Sitter model. It will prove useful to contrast the spheri- 
cal evolution with that predicted by linear theory. We will 
use D(z) to denote the linear density perturbation growth 
factor, normalized so that D(z = 0) = 1. For an Einstein- 
de Sitter Universe, D(z) = 1/(1 + z). Note that this makes 
D oc (t/to) 2 ^ 3 . Similarly, the growth of velocities in linear 
theory is given by 



«un(r) = -^/(S7)A(r,t). 



(A19) 



where /(fi) w J7 6 (Peebles 1980). It is a useful exercise 
to verify that, in its early stages (i.e., small development 
angle) , the spherical evolution model does indeed reproduce 
linear theory. 

Consider the evolution of an initially overdense (or, 
rather, bound) shell. Such a shell will initially expand 
slightly slower than the background, this expansion gradu- 
ally slowing to a complete halt, after which it turns around 
and starts to contract. At turnaround, v(r, t) — 0, so 6 r = 7T, 
and the density is 



1 + A(r,tt») = (3tt/4) 2 



(A20) 



Therefore, at turnaround, the comoving radius of a spherical 
perturbation has shrunk by a factor of (37r/4) 2//3 = 1.771 
from what it was initially. Had the perturbation evolved 
according to linear theory, then turnaround would happen at 
that redshift when the linear theory prediction Au n , reaches 
the value <5 ta : 



A Un (z ta ) = <5ta = (3/5)(3tt/4) 2/3 « 1.062. 



(A21) 



Full collapse is associated with 6 r = 27r. At this time, the 
linearly extrapolated initial overdensity reaches the thresh- 
old value 8 C , 

A lin (z c ) = 5 c = (|) (y)^ 3 « 1-686. (A22) 

This makes it straightforward to determine the collapse red- 
shift z co ii of each bound perturbation directly from a given 
initial density field. In terms of the primordial field linearly 
extrapolated to the present time, A; in>0 , the collapse redshift 
2 co ii may be directly inferred from 



1 + Zcoll = 



(A24) 



1.686 

Formally, at collapse, the comoving radius is vanishingly 
small (7?.(27r) = 0). In reality, the matter in the collapsing 
object will virialize as interactions between matter in the 
shells will exchange energy betwen the shells and ultimately 
an equilibrium distribution will be found. Therefore, it is 
usual to assume that the final size of a collapsed spherical 
object is finite and equal to its virial radius. For a perfect 
tophat density, the object's final size Ra n is then w 5.622 
times smaller than it was initially (Gunn & Gott 1973), i.e. 



Rfin/Ri ,c< 

where Ri 



(18tt 



2x1/3 



5.622 , 



(A25) 



Ri{a co u/ai). 



A7 Underdensities and shell-crossing when f7 = l 

Underdense spherical regions evolve differently than their 
overdense peers. The outward directed peculiar accelera- 
tion is directly proportional to the integrated density deficit 
A(r,t) of the void. In the generic case, the inner shells "feel" 
a stronger deficit, and thus a stronger outward acceleration, 
than the outer shells. 

Once again, to better illustrate our argument, we will 
now specialize to the case of an Einstein de-Sitter model. 
The density deficit evolves as 

9 (sinh6 r - 6 r ) 2 



1 + A{r,z) 



2 (cosh6 r -l) 3 



(A26) 



In comparison, the corresponding linear initial density 
deficit Au n (z): 



A Un (z) = 



Ai,„,o _ _3 2 / 3 (sinh6 r -6 r ) 2/3 
1 + 2 ~ 4 5/3 



(A27) 



The (peculiar) velocity with which the void expands into its 
surroundings is 



. 3 sinh6 r (sinh6 r - 6 r ) 
vM=g ^2 (cosh6 r -l) 2 



(A28) 



As a consequence of the differential outward expansion 
within and around the void, and the accompanying decrease 
of the expansion rate with radius r, shells start to accumu- 
late near the boundary of the void. The density deficit \ A(r) \ 
of the void decreases as a function of radius r, down to a min- 
imum at the center. Shells which were initially close to the 
centre will ultimately catch up with the shells further out- 
side, until they eventually pass them. This marks the event 
of shell crossing. The corresponding gradual increase of den- 
sity will then have turned into an infinitely dense ridge. From 
this moment onward the evolution of the void may be de- 
scribed in terms of a self-similar outward moving shell (Suto 
ct al. 1984; Fillmore & Goldreich 1984; Bertschinger 1985). 
Strictly speaking, this only occurs for voids whose density 
profile is sufficiently steep, since a sufficiently strong differ- 
ential shell acceleration must be generated. This condition 
is satisfied at the step-function density profile near the edge 
of a tophat void. 

For a tophat void in an Einstein-de Sitter Universe the 
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shells initially just outside the void's edge pass through a 
shell crossing stage at a precisely determined value of the 
mass shell's development angle O r = O sc , 



sinhOsc (sinhOsc — B sc ) _ 8 



(cosh © sc - 1) 



so 9 S 



3.53. 



(A29) 



At this shell-crossing stage, the average density within the 
void is 



1 + A(r,t) = 0.1982 



(A30) 



times that the cosmic background density. This means that 
the shell has expanded by a factor of (0.1982)~ 1/3 « 1.7151. 
In comparison, the underdensity estimated using linear the- 
ory at the time of shell crossing is 



(also see Suto et al. 1984; Fillmore & Goldreich 1984). On 
the basis of this, Blumenthal et al. (1992) attempted to re- 
late the derived void characteristics to the observed galaxy 
distribution. Dubinski et al. (1993) (also see Van de Wey- 
gaert & Van Kampen 1993) showed that when this was done, 
then the spherical tophat model provided a rather good de- 
scription of void formation and evolution in their numerical 
simulations. The spherical model is equally succesfull in de- 
scribing the evolution of spherical voids with more generic 
density profiles, and can be employed to demonstrate that 
they will often quickly evolve towards a tophat configura- 
tion (see Fig. 3). Therefore, a description of void evolution 
which is based on the spherical evolution model, a strategy 
which we will follow in the main text, is amply justified. 



A; ln (z sc ) = S v = - ^-^ 



3\ 2/3 (sinh9 sc -O s 



-.2/3 



5/3 



2.81 .(A31) Ag Beyond the spherical model 



In terms of the primordial density field, the shell-crossing 
redshift z sc of a void with (linearly extrapolated) density 
deficit A(jn,o may therefore be directly predicted. For an 
Einstein-de Sitter Universe it is 

1 + Zsc = . (A32) 

2.8059 v ; 

And at shell-crossing, the void has a precisely determined 
excess Hubble expansion rate: 



H s = (4/3) H u (t sc ), 



(A33) 



with H u = H u (t sc ) the global Hubble expansion factor at 

tsc- 

For a spherical underdensity, the instant of shell cross- 
ing marks a dynamical phase transition. It is as significant as 
the full collapse stage reached by an equivalent overdensity. 
Also, as with the collapse of the overdensity the timescales 
on which this happens are intimately related to the initial 
density of the perturbation. The instant of shell crossing 
is determined by the global density parameter fli, the ini- 
tial density deficit Ai of the shell, and the steepness of the 
density profile. In turn, this link between the initial void 
configuration and the void's shell crossing transition epoch 
paves the way towards predicting the nonlinear evolution of 
the cosmic void population on the basis of the primordial 
density field. 

A8 Beyond shell-crossing 

Virtually all early studies of void evolution concentrated on 
analytically tractable configurations of symmetric holes in 
a uniform background, either with or without compensat- 
ing ridges. This allowed Hoffman & Shaham (1982) to argue 
that voids should indeed be a seen as a natural outcome of 
a dissipationless clustering scenario, evolving from deep un- 
derdense regions in the primordial density field. This was fol- 
lowed by a variety of similar numerical studies (Peebles 1982; 
Hoffman, Salpeter & Wasserman 1983; Hausman, Olson & 
Roth 1983). The most extensive and systematic study of 
spherical void evolution, the work by Bertschinger (1983, 
1985) for voids in an Einstein-De Sitter universe, concluded 
that in most viable circumstances voids would develop a 
dense surrounding shell. Following shell crossing at these 
void boundaries, the void would enter a phase of nonlinear 
evolution characterized by a self-similar outward expansion 



We have concentrated on the evolution of spherical pertur- 
bations. However, generic peaks in Gaussian random fields 
are triaxial (BBKS 1986), so it is worth spending a little time 
discussing the evolution of ellipsoidal perturbations. It turns 
out that, as the underdense ellipsoid expands, the spherical 
model becomes an increasingly good approximation. 

A simple approximation for the gravitational potential 
in the immediate vicinity of a density minimum is a sec- 
ond order scheme, which approximates the density field by 
an ellipsoid of uniform density. The evolution of Zoui-density 
regions may therefore be approximated via the equations 
of motion for a homogeneous ellipsoid. The description of 
a void's evolution is therefore analogous to the equivalent 
description of the collapse of overdensities (Icke 1972, 1973; 
White & Silk 1979). Bond & Myers (1996) noted that it 
is possible to incorporate external (anisotropic) influences 
through the appropriate modification of the equation of 
motion (this same scheme was adopted by Eisenstein & 
Loeb 1995). 

In the case of overdense regions, from which collapsed 
halos form, the ensuing nonlinear evolution tends to strongly 
amplify these initial departures from sphericity (Lin, Mestel 
& Shu 1965). The collapse of overdensities typically pro- 
ceeds in an anisotropic fashion, progressing through an in- 
creasingly flattened and elongated configuration before the 
ultimate collapse along all directions is complete. The key 
towards understanding this tendency is the anisotropic force 
field corresponding to nonspherical objects. In the case of 
an overdensity, the effective gravitational forces are directed 
inward, which, in combination with their anisotropy, trans- 
lates into an increased rate of collapse along the shortest 
axis. In the cosmological context, this explains the existence 
of filamentary and sheetlike structures on Megaparsec scales 
(Icke 1973; White & Silk 1979). 

On the basis of the same arguments, voids be- 
come increasingly spherical as they evolve (Icke 1984; 
Bertschinger 1985). That is, the anisotropic peculiar force 
field directed outward will induce the strongest acceleration 
along the shortest axis, causing the void to expand fastest 
along that direction. In contrast, a weaker acceleration along 
the longest axis leads to a smaller rate of excess expansion. 
Hence, the tendency of underdense regions to nullify initial 
asphericities and evolve into an ever more spherical geom- 
etry. Moreover, for a broad range of initial density profiles, 
voids will develop into objects with a distinctly tophat con- 
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figurations. The reason for this is the same as for the spher- 
ical underdensities studied above. This evolution towards a 
tophat profile was indeed observed by Van de Weygaert & 
Van Kampen (1993) for voids in more generic circumstances. 
The homogeneous interior density goes along with a uniform 
velocity divergence. Thus, generic primordial underdensities 
appear to evolve into "super-Hubble expanding bubbles" 
(Icke 1984; Van de Weygaert & Van Kampen 1993). 

Of course the ellipsoidal model has serious limitations. 
It disregards important aspects like the presence of sub- 
structure. More serious is its neglect of any external influ- 
ence, whether secondary infall, "collision" with surrounding 
matter (neighbouring expanding voids!), or the role of non- 
local tidal fields. Yet, it is interesting that in the case of 
voids the homogeneous ellipsoidal model becomes a better 
approximation over an ever increasing volume of space, as 
time proceeds. This has been confirmed by N-body simula- 
tios of void evolution in realistic clustering scenarios which 
show how the matter distribution in the central region of 
(proto)voids flattens out as they expand and get drained 
(e.g. Van de Weygaert & Van Kampen 1993, Fig. 31). On the 
basis of the spherical model (see e.g. Fig. 3) one may readily 
appreciate this generic flattening of the density profile and 
its outwards expansion. Thus, we have the ellipsoidal model 
providing the argument for the sphericity of voids, and the 
spherical model demonstrating why the required conditions 
for the applicability of the ellipsoidal model are generically 
encountered in the case of voids. 



APPENDIX B: RANDOM FIELD 
CHARACTERISTICS 

When evaluating the statistics of a three dimensional ran- 
dom field of density perturbations filtered on a specific spa- 
tial scale R, the spectral moments, 



(Bl) 



play a key role. Here P(k) denotes the power spectrum of the 
unsmoothed density initial fluctuation field, extrapolated us- 
ing linear theory to the present time, and W(kR) represents 
the shape of the filter. For example, if the density field is 
smoothed with a tophat or Gaussian filter, then W(x) is 

_3\ 



Wth(x) = 
W G (x) = 



(3/x ) (sin x — x cos x) . 
exp(-:r 2 /2) , 



(B2) 
(B3) 



respectively. The total volume enclosed by these filters is 
Vth = 4tt_R 3 /3 and V G = (2iv) :i/2 R 3 respectively. The mass 
within the filter is to = pV(l + S), where S represents the 
overdensity: it is this quantity which fluctuates from one 
position to another. If the density fluctuations are small ev- 
erywhere, then the mass within a filter is approximately the 
same everywhere: to « pV. 

In models of hierarchical structure formation, the ini- 
tial fluctuations around the mean density are indeed small. 
Therefore, the correspondence between mass and filter scale 
m oc R 3 suggests that if one wishes to model (proto) objects 
of mass m, one should study the initial density fluctua- 
tion field when it is smoothed on (comoving) spatial scale 
R oc to 1 / 3 , with the exact coefficient depending on filter 
choice. 



Thus, we may analyze any fluctuation field quantity in 
terms of its spatial scale R, or mass scale m. To illustrate, 
consider a power-law power spectrum P(k), 

(B4) 



P(k) oc k n , 
for which 

2/ n -(n+3+2j)/3 
Oj (to) oc to v j " . 



(B5) 



Notice that if — 3 < n < 1, then oo{m) is a decreasing 
function of to. This remains true for any spectrum whose 
"generalized" power spectrum slope, 

dlogP(fc) 



n(k) = 



cZlog k 



(B6) 



is within the range ~3 < n(k) < 1, even if it is not neces- 
sarily constant over the whole spectral range. 

The quantity a (m) quantifies the root mean square 
amplitude of density fluctuations on mass scale to. It will 
feature prominently in this work. In hierarchical scenarios, 
<7o(m) is a monotonously decreasing function of scale, so it 
will serve as a "dimensionless" parameter which character- 
izes and the "scale" of density fluctuations. We will often 
use S(tn) to denote a 2 (m). 

To evaluate the number density of peaks of a scale to 
in the initial Gaussian density field, one must take into ac- 
count the shape of the power spectrum, through the spectral 
parameters i?* and 7, 



R* 



01 



7 



2 

01 



0002 



(B7) 



where ao, 01 and <T2, which depend on the shape of the power 
spectrum, are defined by equation (Bl). 

Bl Simple Peaks Model 

The number density of peaks (minima) of height (depth) 
8 P in an initial Gaussian density fluctuation field, smoothed 
with a filter scale R and corresponding mass scale to = pVf, 
has been been worked out by BBKS. They showed that the 
density of peaks (minima) of height (depth) 



(B8) 



(JQ (to) 



(B9) 



in a Gaussian random field is 

„ (lA Av - H 7 ' c M-»m G(7,7^ 1/2 ) du 
rwWto-yJ— {27T)3/2R s 2 V 

with the spectral parameters R* and 7 given in Eqn. B7. 
For a power-law P(k), P(k) oc k" , some of these integrals 
diverge if one uses a tophat filter. Without loss of physical 
meaning, it is therefore preferrable to use a Gaussian filter. 
The function G is then given by 

-(^-y) 2 /2(l-7 2 )] 



G(7,») 



f 

Jo 



dx f(x) 



exp[ 



VMi-7 2 ) 

in which the function f(x) is defined as 



m = 



(x :i - 3x) 



erf 



5X 1 / 2 



(!) 



+ erf 



(I) 



1/2 



(I) 



1/2 



Six \ 

— + - ) exp(-5a ; 2 /8) 
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-(^--|)exp(- 5a; 2 /2) 



APPENDIX C: FIRST CROSSINGS AND 
LINEAR BARRIERS 

Let f(x) denote the first crossing distribution of a single 
barrier of constant height B: 



Notice that 



B 2 \ dx 



(CI) 



The Laplace Transform of this distribution is 
L(t,B)=J^ dx f(x, B) e~ tx = exp (- V2tB 2 ^) . (C2) 

If we set 5 equal to ao(m) defined by equation (Bl) then 
the distribution f(S,8 c ) gives the excursion set approach's 
approximation for the fraction of mass which is bound up in 
collapsed objects of mass m(S) (Bond et al. 1991). There- 
fore, if n(m) denotes the number density of such collapsed 
halos, then 



p d hi to 



(C3) 



where p is the background density. In models of hierarchical 
clustering, S decreases monotonically with increasing m. If 
the initial spectrum of fluctuations was a power law, then: 
81/ S oc (m/m,) (n+3)/3 with -3 < n < 1. 

By analogy, the fraction of mass in voids which each 
contain mass m(S) is given by 

T(S, 8 V , 8 C ) = f(S, 8 V )- f dsf(s, 8 C , 8 V ) f(S, 8 v \s, 8 C ) :(C4) 
Jo 

the first term on the right hand side is the first crossing 
distribution of the barrier 8 V , and the second term subtracts 
from it the subset of trajectories which had crossed <5 C before 
ever reaching 8 V . Since f(S, 8 v \s, 8 C ) = f(S — s,8 c — S v ), the 
Laplace Transform of T is 



£(t,8 v ,8 c ) = L(t,8 v ) 



ds !F(s, 8 C , 8 V ) e 



f 

J s- 



-t(S-s) 



< I dSf(S-s,8 c -8 v )e~ 

>S-s 

= L(t, S v ) — £(t, 8 C , <5 V ) L(t, 8 C — S v ). (C5) 



Whereas the actual distributions are related by convolu- 
tions (equation C4), the Laplace transforms simply multi- 
ply: in this respect, the Laplace transforms behave similarly 
to one's intuition about the independence of walks to and 
between the barriers. By symmetry, 

C(t,8 c ,8 v ) = L(t,8 c ) - C(t J 8 Vj 8 c ) L(t,8 v - 8c). (C6) 
Inserting this into equation (C5) yields 

L(t,8 v ) - L(t,8 c ) L(t,8 c - 8 V ) 



£(t,S v ,S c ) = 



l-L(t,8 v -8c)L{t,8 c -8 v ) 
sinh[a c x/2l] 



sinh[(<5 c - 8 v )y/2i] ' 



(C7) 



Inverting this Laplace Transform yields equation (1) in the 
main text. 



L(0,8 C ) = 1 



and 



£(0,8 V ,8 C ) = 



(C8) 



8 C — S v 

The first equality shows that all random walks cross 8 C , and 
is interpretted as indicating that all mass is associated with 
gravitationally bound halos. In constrast, the second equal- 
ity shows that only a fraction of all random walks cross 
<S V without first having crossed <5 C ; evidently only a fraction 
/void = 8 C /(S C ~ 8 V ) of the mass is associated with voids. 

Although we do not use this fact in the main text, the 
calculation above can be generalized to include barriers of 
the form B c = 5 c — f3S. These are barriers which are not con- 
stant; rather, their height decreases linearly with S if f3 > 0, 
and increases if f3 is negative. The first crossing distribution 
of linear barriers is Invers e Gaussian; t he associated Laplace 
transform is exp[/3c5 c - ^(2* + (3 2 )} (e.g. Sheth 1998); it 
reduces to equation (C2) when (3 = 0. 

If both barriers change linearly with S, but they have 
the same slope, £> c = S C ~/3S and B v — 8 V ~/3S, then exactly 
the same reasoning which led to equation (C7) yields 



£{t,8 v ,8 c ,/3) =exp(<5 v /?) 



sinh[4V/3 2 +2t] 
sinh[(5 c -S v )y//3 2 + 2t] 



(C9) 



for the Laplace transform of the distribution which 
crosses B v without first crossing B c . (The first cross- 
ing of 8c without crossing 8 V is given by interchanging 
'c' and 'v' in the expression above.) Note that £(0) = 
cxp(<5 v /3) sinh(<5 C| 0)/sinh[(<5 c — 8 v )f3]; as in the constant bar- 
rier model, only a fraction of walks cross one barrier with- 
out first crossing the other. Inverting this Laplace transform 
yields 



SF{S,8 V ,8 C ) = e^]T 



QVP) 2 sin(jTrP) 
v jir 

(jirV) 2 + 2 8 2 



x exp 



where v = 8%/S and V = |<5v|/(5 c - 8 V ). 



2v 



(CIO) 
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